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Abstract 

The objective of this paper is the study of the equilibrium behavior of 
a population on the hierarchical group S7jv consisting of families of indi- 
viduals undergoing critical branching random walk and in addition these 
families also develop according to a critical branching process. Strong 
transience of the random walk guarantees existence of an equilibrium for 
this two-level branching system. In the limit N — > oo (called the hier- 
archical mean field limit), the equilibrium aggregated populations in a 
nested sequence of balls ' of hierarchical radius I converge to a back- 
ward Markov chain on R+. This limiting Markov chain can be explicitly 
represented in terms of a cascade of subordinators which in turn makes 
possible a description of the genealogy of the population. 
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1 Introduction 



Spatial branching processes involve two basic mechanisms, spatial migration 
and branching. These two mechanisms act in opposite directions: the branch- 
ing causes fluctuations of the local population densities which are counteracted 
by the smoothing effect of the migration, and a transient migration is needed 
to sustain an equilibrium of a geographically extended population where each 
individual has an offspring of mean one. Multilevel branching systems (see, e.g. 
DH , GHW], |Wuj 1 involve branching, that is death and replication, at a col- 
lective level. For example, in two-level branching systems both individuals and 
families (that is, collections of individuals that trace back to a common ances- 
tor in the individual branching) reproduce independently. In such systems, the 
fluctuations of the population densities are substantially enhanced compared to 
systems with branching on the individual level only, and in two-level branching 
systems a strongly transient migration is needed to sustain an equilibrium. It 
is well known that Euclidean random walks are transient if and only if the di- 
mension is bigger than 2, and strongly transient if and only if the dimension is 
bigger than 4. In this sense, 2 is the critical dimension for one-level branching 
systems, and 4 is the critical dimension for two-level branching systems. 

In the present paper, we will focus on spatial models with a hierarchical 
(or ultrametric) geographical structure [N islands (blocks of radius one) per 
archipelago (blocks of radius 2), N blocks of radius I per block of radius £ + 1, 
I > 1, cf. SF ). The migration process then is modelled by so called hierar- 
chical random walks: at a certain rate depending on £, an individual jumps to 
a randomly chosen island in distance £. This ultrametric structure leads to a 
separation of time scales as N — > oo, and makes the models particularly suitable 
for a thorough analysis of equilibrium states and cluster formation. It turns out 
that in the hierarchichal mean field limit (with order of N individuals per island 
and N — ► oo) there is a separation of time scales in which the population den- 
sities in the blocks of different radii evolve. For a block of radius £, the natural 
time scale turns out to be N l in the case of one-level branching (see DG2 ) 
and N 1 / 2 in the case of two-level branching. On this time scale, the population 
density in a block of radius £ performs, as N — ► oo, a process whose fluctuations 
are governed by the branching and whose drift is given by a flow of emigration 
and immigration from the surrounding block. For a sequence of nested blocks, 
this leads to a hierarchy of branching equilibria whose structure we describe in 
the next subsection. For the case of two-level branching, the convergence of the 
population densities in nested blocks towards this hierarchy as N — ► oo is stated 
in Theorem 12.4. II and proved in section 

Generically, in our hierarchical model the migration process that sustains 
an equilibrium is at the border to recurrence in the case of one-level branching, 
and at the border to weak transience in the case of two-level branching, as 
N — » oo. In this sense, the hierarchical one-level branching equilibria studied 
in |DG2| correspond to a situation "near dimension 2", and the hierarchical 
two-level branching equilibria studied in our paper correspond to a situation 
"near dimension 4" . Dimension 4 is of considerable interest because it serves as 
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a critical dimension not only for the two-level branching systems studied in this 
paper but also for a number of phenomena in statistical physics including the 
large scale fluctuations of ferromagnetic models at the critical temperature. 

The structures of the family clusters in equilibrium can be best understood 
in terms of the genealogy of the branching system, see |DP 11 IDP21 ISt Wj ) . We 
will describe the genealogy of the equilibrium population in the mean-field limit 
using a cascade of subordinators. 

2 Overview 

2.1 Hierarchies of one— and two-level Feller branching dif- 
fusions 

Consider a large population whose size is fluctuating because of critical re- 
production, and which is subject to emigration of individuals and immigration 
from a surrounding (still larger) reservoir of individuals. The immigration rate is 
given by the population density in the environment, which fluctuates on a slower 
time scale. Now consider an infinite hierarchy of such populations Vg, £ = 1, 2, .., 
where Vg+i acts as an environment for Vg, and think of an equilibrium situation. 
We will study two situations where there is a sequence of time scales such that, 
in the limit of large local population sizes, on the time scale associated with Vg 
the population density Q+\ of Vg+\ remains constant, and given £^ + i = a, the 
dynamics of the population density Q of Vg is of the form 

dQ(t) = dMg(t) - cg((g(t) - a) dt. (2.1.1) 

Here eg is a positive constant which describes the migration intensity into and 
out of Vg, and Mg is a martingale describing the fluctuations of Q. 

In subsection 12 . 2 . 31 we will describe a situation in which the martingale Mg 
has quadratic variation 

d(Mg){t) = Qg{t)dt, (2.1.2) 
hence in this case (|2.1.1|) specializes to 

d(g(t) = VW)dWg(t) - cg((g(t) - a) dt, (2.1.3) 

where W\,Wi, ... are independent Wiener processes. For each i ()2. 1.3(1 is the 
stochastic differential equation of a subcritical Feller branching diffusion with 
immigration f |ANI IEK| ). 

Later on, we will consider a dynamics of the population density Qg which is 
again of the form 1)2.1.1(1 but where the fluctuations are governed by a "family 
structure" of the population. More precisely, the martingale Mg has quadratic 
variation 

d{Mg)(t)=([ x 2 &{t,x)dx) dt. (2.1.4) 

U(0,oo) J 
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where £i(t, x) measures the rescaled number of families of size x within Vi . The 
link between the population density Q(t) and the density £i(t, x) of family sizes 
is given by 



C/(*)= / x£t(t,x)dx. (2.1.5) 
J (0,00) 

The form (|2.1.4|l of the fluctuations of Q indicates that we are dealing with Feller 
branching diffusions of families rather than individuals. This family branch- 
ing shows up in the dynamics which is described by an absolutely continuous 
measure- valued process with density £e(t, x) satisfying the stochastic partial dif- 
ferential equation (SPDE) 

^Zt(t,x) = y/b{t, x) W t (x) + G*&(t, x)-c (-JL.xtt{t, x) + a6' (x)j (2.1.6) 

with c = ci where G* is the adjoint of the generator G of a critical Feller 
branching diffusion given by 

Gf = l^f, (2.1.7) 

W is space-time white noise and 8' is the derivative (in the sense of Schwartz 
distributions) of the (^-function at 0. An equivalent formulation of (|2.1.6|) is 

(2.1.8) 

where 6 e is the Dirac measure in e, and / : (0, 00) — ► R has bounded first and 
second derivatives with f(x) < const x. Note that the first term in (|2.1.(j|l (and 
(|2.1.8fl ) comes from the familiy branching, the second comes from the individual 
branching, and the "mean reversion" term comes from the individual emigration 
and the immigration of infinitesimally small families. 

In addition, (|2.1.6|l shows that the sizes of the single families develop in- 
dependently according to individual subcritical Feller branching diffusions. We 
will therefore call I — 1, 2, . . ., a hierarchy of two-level branching diffusions. 

Two-level branching diffusions have been introduced by Dawson and Hoch- 
berg |DH| as superprocesses over Feller branching diffusions, where all mass 
accumulating in is removed. Therefore, these processes take their values in 
the measures on (0, 00). In fact, it turns out that for t > they have absolutely 
continuous states. In our context, in addition to the set-up of |DH| . there is a 
"continuous immigration of small families". We will see how this fits into the 
framework of immigration processes from the boundary studied in Li and Shiga 
|LS| . For general background on superprocesses and related stochastic partial 
differential equations, see jD] and |DP2| . 

The hierarchies of branching equilibria considered in our paper are motivated 
through a spatial picture which we describe for the case of l|2.1.2l) ("one-level 
branching" ) in subsection 12.2.31 and for the case of (|2. 1.411 ( "two- level branch- 
ing" ) in subsection 12.41 The case of a hierarchy of one- level branching systems 
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was studied by Dawson and Greven |DGHlUU2"] in the context of super-random 
walks (or interacting Feller branching diffusions). 

For any given 9 > (which in the geographical model will play the role of a 
"global population density" ) we will construct the hierarchy 

(2.1.9) 

in terms of an entrance law for a backward Markov chain where the conditional 
law of Q given Ce+i = a is an equilibrium state of l|2.1.1|) . More precisely, in 
subsection 14 . 1 1 we will show the following result. 

Proposition 2.1.1 Let (q)^>i be a sequence of positive numbers, and let us 
distinguish two cases: 

a) Assume ^gcj 1 < oo. For a > 0, £ E N, let H e (a, .) be the equilibrium 
distribution of 

b) Assume cj 2 < oo. For a > 0, £ E N and c = q, let £^ i0 be an equilib- 

_______ ( 2^ 

ram state of \2.1.b}) . and II^ (a, .) be the distribution of J^ Q ^ x^i ta {x)dx. 

In both cases, for every 9 > the backward Markov chain with transition 
probability 

P{QeA\Ci +1 =a)=Il i (a,A), 
where Id is either orY\^ 2 \ has a unique entrance law {Qf }e=....2,i satisfying 

EC- = 9 (2.1.10) 

and 

lim d = a.s. (2.1.11) 
2.1 A A cascade of subordinators 

To work out parallels between the one- and two-level branching situations de- 
scribed in subsections 12.2.31 and 12.41 and to discuss aspects relevant for the 
genealogy of the hierarchical branching equilibria, we write He(a, .) for the equi- 
librium distribution of l|2.1.1|) in the two cases (|2.1.2|l and <|2.1.4|) (which corre- 
spond to cases a) and b) in Proposition 12.1.1(1 . 

In both cases the parameter a enters as a factor into the immigration rate of 
a continuous state branching process, hence Hi(ai + 02, .) = n^(ai, .) *Ibj(a2, .). 
Therefore the Tle(a, .) are infinitely divisible distributions on (0, 00) and there 
exist subordinators (that is, processes with stationary independent non- negative 
increments) Si(a) , a > 0, such that 

C(S e (a)) = U e (a, .). (2.1.1) 

We will see in subsection 16 . 21 that in case a) the Si are Gamma processes. 

In both cases, the Levy-Khinchin decomposition (see |K2| . Chapt. 15) 
S((a) — ^j^Hi describes (asymptotically as N — > 00) the partitioning of the 

equilibrium population in B^ N ^ into subpopulations stemming from one and 
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the same immigrant into , given that the population density of the sur- 
rounding block is a. 

Recall from ((2.1.1(1 and Proposition ^. 1 . ll that 

xJl e (a,dx) = a. (2.1.2) 

Therefore, 

ES e (a)=a. (2.1.3) 

Let us denote (in either of the two cases) the Levy measure of Si by pn, 
and the second moment of hi by mi. An iteration of the Levy-Khinchin repre- 
sentation (which can be interpreted in terms of the genealogy of the branching 
hierarchy, see subsection 14.3(1 will show that the Levy measure of the iterated 
subordinator Si(Si+±(. . . Sj—±)) has second moment mi + . . . + mj—i (see sub- 
section Under the condition 

oo 

Y^m^ < oo. (2-1.4) 

1=1 

we will prove in subsection 14 . 21 that for each 9 > the limit in distribution 

Q 9 t = d- lim Si(Si +1 . . . (2.1.5) 

exists, has expectation 9 and defines an entrance law with respect to (IL(a, .)). 
In particular one has 

$= d S t (S t+1 ...S j - 1 {$))), j>L (2.1.6) 

For each j > I this gives a decomposition of Q, which asymptotically as N — > oo, 
stands for the partitioning of the equilibrium population (g N ' 8 ^ in Bg N ^ into 

subpopulations stemming from one and the same immigrant into B^\. 

The summability condition (|2.1.4|l is equivalent to the transience condition 
(|2.2.4|l in the one-level case, and to the strong transience condition (|2.3.10|l in 
the two- level case, since we will show that mi — 1/2q in the one-level case, and 
mi = l/4c| in the two- level case (see Remarks 16.2.31 and 13.2.7(1 . 



2.1.2 Genealogy 

In Section l4.3l we develop a genealogy of the jumps occurring in the cascade of 
subordinators. The idea is that given a jump of Se+i(-) at time ti there will be 
a set of jumps of Si(-) that occur in the time interval {Si+\{ti— ), Se+i(ti)) and 
these level £ jumps will be said to be descendants of the level (£ + 1) jump. In 
subsections ((4.3(1 . ((4.41) and ((4.611 we use this idea to develop the full genealogical 
structure of the population associated with the entrance law. This leads to a 
decomposition of the population into a countable collection of subpopulations 
of individuals having a common ancestor. For the case of critical individual 
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branching this was done in DG2 , for the two-level branching case this is new. 
We will work out the parallels between the two cases in a general framework, 
which also sheds some new light on the results of DG2 . 

Intuitively this genealogy describes the limiting genealogical structure of 
the spatial branching equilibria with hierarchical geographic structure described 
in the Introduction as the parameter N — > oo and the analogue of the "clan 
decomposition" of the equilibrium of super-Brownian motion (e.g. [DPI ). 

2.2 Hierarchical geography, random walks and branching 
equilibria 

2.2.1 A class of random walks 

In order to give a precise formulation for the spatial system we now describe 
the set iljy of sites on which the spatial population lives. For fixed N £ N, 
let f2jv be the countably many leaves of a tree all of whose inner nodes have 
degree N + 1. In other words, each node in depth £ + 1, £ = 0, 1,..., has 
N successors in depth £ and one predecessor in depth I + 2. For two sites 
y, z G fijv, their hierarchical distance dpf(y,z) is defined as the depth of their 
closest common ancestor in the tree. Note that djv is an ultrametric, that is, 
djv(y, z) < maji{d,N(y, x), cLn{x, z)}. We define the individual migration on flpj 
as follows. Let q[ , q% > ■ ■ ■ be positive numbers with J2e Qe^ < 00 • At rate 
<jfj the individual makes a jump of distance £ (i.e., it waits for an exponentially 
distributed time with parameter y\ qi N ^ and then jumps a distance j with 

probability qf 1 / J2e It ' )■> choosing its arrival site uniformly among all sites at 
a distance £ from its previous site. 

The set ft^ can be identified with the set of sequences in {0, 1, . . . , N — 1} 
almost all of whose entries are equal to zero. With component-wise addition 
mod N, fijy is a countable Abelian group (the so called hierarchical group of 
order N). Note that d^{y,z) is translation invariant; it will be written as 
\y — x\. The migration specified above is a (continuous time) random walk on 
iljv called hierarchical random walk. Hierarchical random walks, in particular 
their transience and recurrence properties, are studied in |DGW2j . jDGWSj . 

2.2.2 A system of branching random walks 

We now introduce a system of branching random walks on £1^- This is given 
by a system of particles undergoing symmetric random walks with migration 
rates q^ N ' together with branching. Branching and migration are assumed in- 
dependent. We specify the branching mechanism as simply as possible: after a 
unit exponential time, an individual, independently of all the others, either is 
removed or splits into two, each case with probability 1/2. 

Remark 2.2.1 ]LMW}/ . ]Gr}j Assume that the migration rates q[ N ^ are such 
that the individual random walk is transient. Then, for each 9 > 0, there exists 
a unique branching equilibrium with a mean number of 9 individuals per site. 
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This equilibrium is of Poisson type, i.e. the equilibrium population $ is the 
superposition $ = Y^i^i °f a Poisson system X^<5*> of families (each family 
consists of a collection of individuals which in a historical picture can be traced 
back to a common ancestor). 



2.2.3 The hierarchical mean-field limit of a branching equilibrium 
"near dimension two" 

Now consider, for a large N, the total number Z^ N \t) of individuals in a closed 
ball of radius I at time t. Note that has N e elements, and look at the 
time evolution of the population density (or block mean) Ai N \t) = zj, N \t)N~ e 
in Bg at time scale tN l . This corresponds to the classical Feller branching 



diffusion approximation EK . However, in order for the immigration into and 
emigration from B^^ to produce a nondegenerate drift term in the limit N — > oo 
one must adjust the migration rates. The appropriate adjustment is 



q M = c t - 1 N-V-», (2.2.1) 
where the eg do not depend on N and satisfy 

lim sup ^±1 < oo. (2.2.2) 

t ct 

We will call a random walk on fl^ with jump rates q\ N ^ of the form H2.2.1|) a 
(1, (eg), N)-random walk. The following proposition is proved in DGW2 . 

Proposition 2.2.2 Consider a (1, (ci), N)-random walk on Qn with rates given 
by \2.2.1]l . Assume that 

lim sup— <N. (2.2.3) 



Then the random walk is transient if and only if 

zZ c e 



^c^Koo. (2.2.4) 



Now assume that (cj) satisfies (|2.2.2|) . Then for each N > limsup^ 



ce+i 
a ' 

because of Proposition 12.2.21 and Remark 12.2.11 there is an equilibrium for the 
system of (1, (q), 7V)-branching random walks with mean 9 for each 8 > . We 
now consider the corresponding equilibrium population densities j4« in B^ N ^ , 

£ = 1,2,... (where here and below (Bi N ') denotes a sequence of nested balls 
of radii t in Qn)- In order to identify the limit as N — » oo of this sequence of 
population densities we must consider the dynamics of A^ in its natural time 
scale N e . Let us first discuss on a heuristic level why, in the limit TV — > oo, the 
drift term on the r.h.s. of (|2.1.1|l arises on the time scale N l for the population 
density in B^ N \ Because of the ultrametric structure of f2jv, an individual in 
B^ has to make a jump of size > £ + 1 in order to leave B^ N \ Because of 
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(|2.2.1[l and (|2.2.2[) . for TV large, jumps of size > £ + 1 happen rarely compared 
to jumps of size £ + 1 (since qi+k = o(qi+i) as TV — * oo for k > 1). Hence the 
individual emigration rate from B) N ' on time scale N e is ci (asymptotically as 
TV — > oo). Concerning immigration into Be N \ again because of (|2.2.1|l . it is 
only the environment in B^fl that is relevant for large TV. An individual in 
B^\ \ B^ N ^ necessarily has to jump a distance £ + 1 in order to make it into 
b\ N \ and on average every (TV — l)-st of these jumps will take the individual 
into b\ N ^ (note that sj+i is TV times as large as £?j )■ Since the block mean 
A^[\ does not change on the time scale N e as TV — ► oo, the total immigration 
rate into B^ on this time scale is (asymptotically as TV — ► oo) of the order 

z£!N*q£l/N = (JV'+^ffi)^ = QAffiiV*. (2.2.5) 

This suggests that the limiting dynamics of the population densities , in 
their natural time scales, as TV — * oo is given by lf2~TTT|) with a = A^\. The 
separation of time scales on balls of different radii that underlies the previous 
discussion is a feature of the hierarchical random walks, which is due to the 
ultrametric structure of (see DGW2 ) . This is also explained in more detail 
in Remark 12.3.41 below. 

Instead of branching particle systems, Dawson and Greven |DG2 consider 
super-random walks (or so-called interacting Feller diffusions) on Note 
also that the definition of the random walk in |UG2| is slightly different but 
asymptotically equivalent as N — * oo to the one used in this paper. In |DG2| 
the sites to which a jump is made are chosen with uniform distribution on a ball 
rather than on a sphere. However, the "interior" of the ball is asymptotically 
negligible as compared to the sphere as TV goes to infinity. 

A particle system analogue of Theorem 4(b) of |DG2j is the following, which 
we state without proof. 

Proposition 2.2.3 Consider a sequence (eg) satisfiying conditions \2.2.J$ and 
\2.2.2)) for transience. For N large enough such that \2.2.ty is met, let the indi- 
vidual migration process be a (1, (q), TV) -random walk. For 8 > let Z^ N,e \t) 
denote the total number of individuals in at time t and A^ N ' \t) :— 

Zg N ' 6 \t)N~ e ' be the population density at time t in B^ in the Poisson type 
branching equilibrium population on Q,n with mean number 9 of individuals per 
site (see Remark \2.2.1\) . Let {Q } be the entrance law provided by Proposition 
\2.1.1\ case a). Then 

{Af'^iO)}^ =>■ {QVn as TV^oo, 

where ==>■ denotes weak convergence of finite dimensional distributions. 

Let us now explain in which sense transient (1, (q), TV)-random walks can 
be interpreted as random walks "near dimension 2" . 
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Definition 2.2.4 'DGW21 Let Z be an irreducible transient random walk on a 
countable Abelian group T. Its degree of transience is defined by 

7 := sup{r? > : E L r ' < oo}, (2.2.6) 

where L is the last exit time of Z from G T. 

Expressed in more analytic terms, the degree of transience of Z is 

/oo 
t v p t (0,0)dt < oo}, 

where pt is the transition probability of Z DGW2] ISaWj . 

For simple symmetric random walk on the d-dimensional lattice Z d , it is well- 
known that dimension 2 is the borderline for transience. For d > 2, the degree 
of transience is d/2 — 1, since the rate of decay of the transition probability is 
Pt(0,0) - const.^/ 2 . 

Remark 2.2.5 \DGW2j a) Let < c < N. Then the (1, (c e ), N) -random walk 
on f2jv * s transient if and only if c > 1 . In this case its degree of transience is 
logc/(log N — logc). Thus for fixed c the transient (1, (cr), N) -random walks on 
riisr have degrees of transience 0(1/ log TV) and therefore asymptotically as N — > 
oo, can be viewed as analogues of random walks "near (Euclidean) dimension 
2". 

b) Assume that (ce) satisfies conditions \2.2.4\ and §2.2.2\ . and put c := 
limsupc£ + i/c£ > 1. Then, for all N > c the (1, (q), N)-random walk on f2jv is 
transient with degree of transience in the interval [0,logc/(log N — logc)]. 

Since certain properties of systems of branching random walks such as per- 
sistence and structure of occupation time fluctuations depend only on the degree 
of transience of the random walks, branching populations whose migration pro- 
cess is a hierarchical random walk can give insight into the behavior of a larger 
class of branching populations whose random walks have the same degree of 
transience. 

2.3 Two-level branching systems in a hierarchical geogra- 
phy 

2.3.1 Strongly transient migration 

Whereas the situation described in subsection l2.2.3l gives an analogue to a sit- 
uation "near dimension 2" , our main focus later on will be on the analogue to 
a situation "near dimension 4" . In this context we will consider the (stronger) 
mass fluctuations induced by a critical reproduction of whole families (of mu- 
tually related individuals), together with a (stronger) smoothing caused by a 
strongly transient migration. 
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Definition 2.3.1 An irreducible random walk Z on a countable Abelian group 
r is called strongly transient if 

E L < oo (2.3.7) 

where L denotes the last exit time of Z from 0. A transient random walk with 
EqL = oo is called weakly transient 

Note that strong transience is equivalent to 

tp t (Q,0)dt < oo 



and a necessary condition is that the degree of transience be equal to or greater 
than 1 DGW2 . Moreover, as mentioned above simple symmetric d-dimensional 
random walk has degree of transience d/2 — 1, and it is strongly transient iff 
d > 4. 

In order to introduce a family of strongly transient random walks on fijv we 
replace ()2.2.1[) by 



IN) 



ct-iN-W 2 , (2.3.. 



where the q do not depend on i and satisfy l|2.2.2|) . We will call a random walk 
on Qjy with jump rates 3/ of the form (|2.3.8|l a (2, (q), N)-random walk. 

Proposition 2.3.2 'DGW2J Consider a (2, (q), N) -random walk on with 
rates given by 12.3. <S)) . Assume that 

l imsup f£±i <N 1 / 2 . (2.3.9) 

i ct 

Then the random walk is strongly transient if and only if 

^2cJ 2 <oo. (2.3.10) 

i 

Remark 2.3.3 WGW2f a) Let < c < N 1 / 2 . Then the (2,(c e ),N) -random 
walk on £In is strongly transient if and only if c > 1 . In this case its degree of 
transience is 

log N + 2 log c 
logTV - 21ogc' 

Thus for fixed c the strongly transient (2, (c ), N) -random walks have degree of 
transience l + 0(l/logiV) and therefore asymptotically as N — » oo can be viewed 
as analogues of random walks "near (Euclidean) dimension 4"- 

b) Assume that (q) satisfies conditions \2.3.1U$ and \2.2.2fy . and put c := 
limsupc£ + i/cf > 1. Then, for all N > c 2 the (2, (cg),N) -random walk on 
CIn is strongly transient with degree of transience in the interval [1, (logiV + 
21ogc)/(logiV - 2 log c)]. 
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Remark 2.3.4 The natural time scale for the strongly transient (2, (c e ), N)- 
random walk in B^ N ^ is N l l 2 . More 'precisely, as shown in ]DGW2j . asymptoti- 
cally as N — => oo for this random walk on the time scale iV^/ 2 only the migrations 
within the ball B^ and the surrounding ball B^J are relevant. Similarly, this 
occurs for the transient (1, (c ), N)-random walk in time scale N l (see 'DGW21 
for details). This effect is basic for the limiting hierarchy of branching equilibria 
obtained in this paper. 

2.3.2 Two-level branching equilibria 

A main objective of this paper is to study two-level branching systems for a 
migration which is on the border between strong and weak transience - recall 
that a strongly transient migration is required for the existence of a branching 
equilibrium. Thus, for Euclidean random walks, d = 4 is the critical dimension 
for a two-level branching system in the same way as dimension d — 2 is the 
critical dimension for a one- level branching system. 

We are going to study two-level branching systems on Sljy. Consider a sys- 
tem of (2, (q), iV)-random walks on Qn such that {ce) satisfies the conditions 
(|2.3.9[1 and (|2.3.1U|) for strong transience, and recall from Remark 12.2.51 that 
these random walks are close to the border to weak transience at least for 
large N if (a) does not grow superexponentially. Introduce, in addition to the 
individual branching and migration, a family branching: independently of ev- 
erything else, after a unit exponential time each family $i either vanishes or 
reproduces resulting in two identical copies each case with probability 

1/2. After a reproduction event 4>£ and evolve as independent one level 
branching systems. This creates the basic two-level branching system *$>( N '{t) 
which is started with the family system at time to given by W N '(to) — J2i 
described in Remark 12.2.11 

The following result is the analogue of GHW for two-level branching sys- 
tems on fijv- 

Proposition 2.3.5 Assume that 

(i) the random walk on is strongly transient and 

(ii) ty( N >(to) = where corresponds to the family decomposition of 
an equilibrium state for the one-level branching random walk with mean number 
8 of individuals per site. 

Then as to — * —oo, the two-level branching system ^^(O) converges in distri- 
bution to a translation invariant equilibrium ^>^ N ' 9 \0) with a mean number 8 of 
individuals per site. 

Remark 2.3.6 The notation ^( N \t) and ty( N > e \t) will be used throughout to 
denote the two level branching system and the equilibrium process with mean 8, 
respectively. 

Remark 2.3.7 Greven and Hochberg 'GHl have obtained more general condi- 
tions under which the convergence to equilibrium as in Provosition \2.3.o^ occurs 
as well as conditions under which it fails. 



14 



Now we consider a system of (2,(ce),N) random walks with (ct) satisfy- 
ing conditions l|2.3.10|) and l|2.2.2l) for strong transience. Then for each N > 

^limsup^ because of Propositions 12.3.21 and 12.3.51 there is a two-level 

branching equilibrium with a mean number of 9 individuals per site 

for each 9 > 0. A main objective of this paper is to study the equilibrium 
structure that arises in the limit as N — > oo of the corresponding sequence of 
family structures in the blocks -E?| . 

2.4 The hierarchical mean-field limit of a two-level branch- 
ing equilibrium "near dimension four" 

2.4.1 Local normalized family-size process 

Let, for fixed N, ^ N ' e '(t),t £ M, be the equilibrium process of the two-level 
branching system as provided by Proposition 12.3.51 Denote the number of 
families in ^ N ' 6 ^(t) having j individuals in -B« by ni N (t,j) and write 
Hg N (t) := nf* (t, j)Sj for the local family- size process. For each fixed £, 
H^ N ' 6 ^ (0) is a random measure on N which describes the population of equilib- 
rium family sizes in the block B^\ We note that the process {H^ N ' 6 \t)} t i=s. 
can be viewed as a branching Markov chain on Z + with instantaneous killing at 
0, the Markov chain on Z + being a standard subcritical binary branching pro- 
cess with immigration. Note that J2j 3 n t' S \ti j) — (t), the number of 

individuals in s£ at time t . Now consider the equilibrium normalized family 
size process defined by 

4 N > e \t) =Y,N- l ^ e \tN*'\3)5 jN - t/ „ (2.4.11) 
j 

In other words, for < a < b < oo, 

4 N ' 6 \t)(a, b) = N- l ' 2 m {N l 'X N e / 2 ( ai b)). (2.4.12) 

Note that the natural time scale in which to observe the subpopulation in B^ N ^ 
in this case is N^/ 2 and not N e as was the case for one-level branching (see 
Remark EO) . 

For each i, N and t, r{f ,e) (t) is a random measure on (0, oo). More precisely, 
we take as state space the set M 1 (0, oo) of Radon measures /i on (0, oo) that 
satisfy the condition j xn(dx) < oo. (Note that we do not keep track of families 
of size 0.) 

The corresponding normalized population mass in B{ N ^ (the "radius £ block 
average" ) is given by 

C f '*>(*) = / x4 N > e \t,dx) =Y,N-^ e \NV%j) = N-% N > 6 \NV% 
J j 

(2.4.13) 
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and in terms of ^^'^(t): 




(2.4.14) 



2.4.2 Convergence theorem 

We now state our main result that makes precise the sense in which the entrance 
law described in section 12.11 approximates the two level spatial equilibrium in 
fijv obtained in Proposition ^ . 3 . 5l when the parameter N — ► oo and the random 
walk satisfies (|2.3.1U|) . 

Theorem 2.4.1 (Hierarchical mean-field limit) Consider a sequence (q) satis- 
fying conditions $2.2.2}) and $2.3. 1U\) for strong transience of the (2, (ce), N)- 
random walk. For fixed N G N obeying $2.3.9}) . and a sequence of nested 
blocks m S]/v, let {Q h r ' 6 \o)} i<=n be the radius t block averages (defined in 

\2.4.1?j\j ) of an equilibrium two-level branching system with an expected number 
9 of particles per site. Let {<^f}e^ be the entrance law provided by Proposition 
HTH1 case b). Then 

The proof of this result is based on the spatial ergodic theorem for the 
equilibrium random field on f2jv obtained in section 15.11 a separation of time 
scales property derived in section l5"2l and a diffusion limit theorem for the family 
size processes {r]g N ' 9 \t)} as N — + oo obtained in sect ion IS31 Using these results 
the proof of Theorem l2.4.1l is given in section I5T41 

3 Super subcritical Feller branching 

In this section we continue the investigation of diffusion limits of two-level 
branching populations without geographical structure, which were introduced in 
|DH| . In our case, these are superprocesses whose basic process is a subcritical 
Feller branching diffusion killed at (this killing corresponds to the removal of 
void families) . With a view towards the application to the hierarchically struc- 
tured geographical model, we will concentrate in subsection 13.11 on an initial 
condition of many small families, which in the diffusion limit corresponds (on 
a heuristic level) to an intial condition ooSq. In subsection 13 . 2 . ll we investigate 
time stationary super subcritical Feller branching processes which arise as dif- 
fusion limits of two-level branching populations with a high-rate immigration 
of individuals. The simplest situation is to think of each immigrant individual 
founding a new family; in the diffusion limit this leads to super subcritical Feller 
branching diffusions with immigration of oo<5o at a constant rate (abbreviated 
by SSFBI). Again with a view towards the geographical model, we will consider 
the situation where (only) every once in a while a newly immigrant individual 
belongs to an already existing family. If this happens relatively rarely, then the 
diffusion limt remains to be SSFBI, see Proposition 13 . 2 . ll and Corollary 13. 2. 21 
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3.1 Diffusion limit of two-level branching particle systems 

Fore > Oande > 0, consider the M 1 (0, oo)-valued family-size process {H £ (t, dx)} 
of a two-level branching particle system (without geographical structure) with 
branching rates equal to 1/e at both levels and subcritical at the individual level 
with subcriticality parameter sc. (An example is the local family size process 
ff( N >®) (defined in subsection l2.4.l|) run at time scale N e ^ 2 and with immigration 
suppressed; here, c — ci and e = N~ £ / 2 .) 
Consider the rescaled family-size process 

rf(t, (x u x 2 )) := sH £ (t, { Xl /s,x 2 /e)), t > 0. (3.1.1) 

Proposition 3.1.1 Let rj £ (t) be as in ^J.l.l]) . Assume that rj e (0) — eLfJ^L-J 
where a > and x > are fixed. Then 
(a) as e — > 0, 

{?f (*)}t>o => {£(*)}*>o (3-1.2) 

in the sense of weak convergence of My([0, oo))-valued eddlag processes, and £(t) 
is the Mf (JO, oo)) -valued superprocess starting in £(0) = aS x , whose motion is 
the subcritical Feller branching process with generator G c given by 

1 d 2 d 

2 dx 2 dx 



G cf=\7;x—- cx —)f> (3- 1 - 3 ) 



acting on functions 

f £ C 2 ([0,oo)) := {/ 6 C 2 ([0,oo)) : lim f(x) = 0, /(0) = 0}. 

X — >oc 

(Here M/([0, oo)) denotes the space of finite measures on [0,oo).) 

(b ) The law of the process £ (t) is uniquely determined by the Laplace functional 

as follows: 



E a5m ^cxp(- f(y)at, dy))j = exp(- J V t f(y) £(0, dy)) 

= exp(— au(t, x)), (3.1.4) 
where Vtf(x) = u(t,x) is the unique solution of the non-linear p. d.e. 

^hA = G Cl u(t,x)- l -u(t,x) 2 (3.1.5) 
u(0,x)=f(x). 

Proof This is essentially Theorem 4.1 of DH . 



Remark As an application of (|3.1.4|l and 1|3.1.5|) we obtain the compact support 
property of £ appearing in Proposition 13.1.11 Assume that £(0) has compact 
support and let Rt denote the range of £(s) up to time t. Then following the 
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method of Iscoe as in Theorem 1.8 of |LS| or Theorem A of D LM| one can show 
that Rt is bounded almost surely. This involves showing (as in DLM ) that for 
c > the equation 

1 d 2 u du 1 2 
2 X d^ = XC ~d~x + 2 U ' 

»<»>=»• % 

has for any a > a blow-up at some finite x. 



x=0 



3.1.1 Evolution equation and entrance law 

In Proposition 13.2.11 we will prove an extension of Proposition 13.1.11 which in- 
cludes immigration. In this subsection we obtain some properties of the solution 
of the evolution equation that will be used there. 
Let 

Ci(0,oo) := {/ e C(0,oo), lim f{x) = 0, |/(x)| < const • x} (3.1.6) 

X — >oo 

and (Tt) be the semigroup of the Feller branching diffusion with subcriticality 
parameter c, absorbed at zero. 
Let Vtf be the solution of 



Lemma 3.1.2 Let f £ Ci(0,oo), / > 0. Then for t > 0, V t f(x) is 
tiable at zero and 

(Vjy(0) = J f(y)K t (dy) ~\f j (VJ) 2 (y) Kt _ s (dy) ds (3.1.8) 

where n t is the (Tt)- entrance law given by $6.2.16)) in the Appendix. 

Proof. First note if / 6 0^(0,00) then V t f(x) < T t f(x) so that V t maps 
Ci(Q, oo ) into itself. Using the evolution form of l|3.1.7|l . 

lim *M = lim T -lM _ I f lim T J^im ds (3.1.9) 

elO £ e|0 e 2 J Q elO 8 

The result then follows from (|6.2.17|l . ■ 

Proposition 3.1.3 Let £, £ (i), t > denote the super subcritical Feller branching 
diffusion (without immigration) process starting in e~ 1 S e at time 0. Then as e — > 
0, £ e converges in the sense of weak convergence of M (0,oo) -valued continuous 
processes on the time interval [to, oo) for all to > to a measure-valued diffusion 
£° where, for all t > 0, £°(t) is an infinitely divisible random measure with 
Laplace functional given by 

Eexp(-(£°(t),/)) =exp(-(F t /)'(0)). (3.1.10) 
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Proof. Because Vtf(Q) = we have 

Eexp(-<£°(i),/}) = hmEexp(-(e £ (i),/)) (3.1.11) 

= limE e -i,,exp(-(£(t),/>) 

ej,0 

= hmexp(- £ - 1 (F i /)( £ )) 

= hmexp (s-\(V t f)(e) - (V t /)(0)) 

= exp(-(T4/)'(0)). 

■ 

Remark 3.1.4 Smce £ (t) is infinitely divisible, its Laplace transform must be 
of the form 

(-(?(t)J))=ew(-f(l-e-( m 'fi)K t (dm)) , / e d((0, oo)). 

(3.1.12) 

for some uniquely determined measure K t on M{Q,oo), the space of Radon 
measures on (0, oo). The measure K t is the canonical measure of£°(t). 

We can now put these results into the framework of |LS| . 

A crucial property of the entrance law (k 4 ) given by l|6.2.15[l . which follows 
immediately by partial integration from the density (|6.2.16|) . is given by the 
following lemma. 

Lemma 3.1.5 For all bounded continuously differentiable functions g on [0, oo) 
with 0(0) = 0, 

lim / g{x)K t (x)dx = g'(0). (3.1.13) 



j exp I 



We fix a strictly positive function p G D(G C ) with 

p(x) — x for x 6 (0, — ], p{x) = 1 for x > 1. (3.1.14) 

Note that such a p meets condition [A] in |LS| . We take as state space M p := 
{p G M(0, oo) : / p{x)p{dx) < oo}. 
Following 1LS) we put 

C p (0, oo) := {/ G C(0, oo) : |/| < constp, lim f(x) = 0}, (3.1.15) 

X — *oc 

D p (G c ) := {/ G D(G C ) : /, GJ G C p (0, oo)} (3.1.16) 

and 

/>OG 

K o+(ff):=lim/ g(x)nt(x)dx, g G D P (G C ). (3.1.17) 
Combining i3.1.1()t . (13.1.1^ . (I3.1.13ll and l!3.1.17ll we obtain 

(l-e-< m ^)K t (dm) = (V t f)'(0) = K 0+ (V t f), f G d((0, oo)), (3.1.18) 

where Vt/ is the solution of (|3.1.7fl . 
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3.2 Super subcritical Feller branching diffusion with indi- 
vidual immigration 

3.2.1 Diffusion limit with immigration 

We now extend Proposition ^ . 1 . Ti to include immigration, taking a fixed to as ori- 
gin of time. Since in our application the population from which the immigrants 
come is structured into families that undergo family branching we incorporate 
multitype immigration and label the set of possible families of immigrants by 
7:= [0,1]. 

Let Af 1 (7 x (0, oo)) denote the set of Radon measures /jon Jx (0, oo) satisfy- 
ing J lx , ^ xfi(dy, dx) < oo. We denote the single atom measure corresponding 
to one individual of type y| £ I by 6y* 1 1 . 

Consider the M 1 (I x (0, oo))-valued family-size process {Hf(t,dy,dx)} to <t 
with branching rates equal to 1/e at both levels, critical at the family level 
and subcritical at the individual level with subcriticality parameter ec (i.e. a 
mean offspring number of 1 — ec per branching event) , and with immigration of 
individuals of type y\ S /, given by fc, i, k € N, at rate ca £ k /e 2 with J2k a % = a 
and lim e ^o sup fe a| = 0. (The motivation for this comes from our geographical 
model, we will see in Section [S] that this setting corresponds to the situation 
where the surrounding population, which serves as the source of immigration, 
is thought to have a frozen family structure.) Consider the rescaled process 
{^/Wlto^t defined by 

Vj(t,{y s k } x (x u x 2 )) :=EHf(t,{y%} x A^)). (3.2.1) 

Proposition 3.2.1 Assume that as e — > 0, r)j(t ) =$> ^ £ M 1 (J x (0, oo)) and 
a e — > a in the sense of weak convergence of finite measures on I , where a is 
a nonatomic measure whose total mass a = ot(I) plays the role of the overall 
immigration rate. Then as e — » 0, 

fof(*, •)}*„<* =>{&(*)}*,<« 

in the sense of weak convergence of Af 1 (7 x (0, oo)) -valued cadlag processes 
on the time interval [to,oo), where {£/(£)} is the measure-valued diffusion with 
generator 

<&F(n) = 0})(/i, G c 4>) + </>))(?, 2 } (3.2.2) 

+f'(M))cf$£(v,x)\ ai=0 a(dy), 

G c is the operator given by ys.l.'d}) and D ((£>) denotes the class of functions of the 
form Ff i( f,(p) = /((/x, <f>)) where f s C^(R) and <f> is continuous function on I x 

(0, oo) with d bounded and continuous on I x [0, oo) and | supj, <f>(y, x)\ < 

const • x. 
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Proof. r/ £ is given by a pregenerator 25 £ acting on the class of bounded contin- 
uous functions F on A/ 1 (7 x (0, oo)) given by 



<5 £ F(/i) = J2 -£[Ffr + ety.*) ~ FM] (3-2.3) 

k 

k 3=1 



k 3 = 1 

\{l + ec) EE^(M ~ ^ + **«£,&■-!)*)) - ' J:) 

k 3 = 1 



Here the first term comes from the immigration, the second and third from the 
critical branching at the family level and the fourth and fifth from the subcritical 
branching at the individual level with subcriticality parameter c > 0. 
For F S D(&), <5 e F takes the form 

© e F(X) (3.2.4) 

= E^[/(^>>+ e ^ e ))-/(^>)] 

k £ 



2 

fc i=i 

1 ^ 

+ 2 (1 ~ ec) EE^' ^> - ^(W.ie) + (3 + l)e))) - /((m, <t>))\ 

k 3 = 1 

e 2 

1 

+ 2 (1 + ec) E £[/«/*, 0) - z<t>(Vk,je) + s<t>(Vk, (j - IK))) - /((M, <£»] 

3 = 1 

e 2 

Tightness of the family {i]j}o <e <i is proved by a standard argument as 
in |Wu| or |DZj . Using a Taylor expansion for the functions / and <f> it can 
be veri fied th at as e -> 0, for F/,^ G D(0), e F/^(/x) -> <5F W . Then by 
Lemma lb . 2 .4l below. we obtain bounds on the third moments uniform in s. Using 
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this and the tightness it follows that any limit point of the laws of the family 
{r]j} satisfies the martingale problem associated to the generator 25. Finally a 
standard argument (e.g. |DP2| . proof of Theorem 1.1) shows that any solution 
of this martingale problem has Laplace functional given by 

E Mo exp(- f(u,x)£i(t,dy,dx)) (3.2.5) 



= exp(- J V t f(y,x)i4)(dy,dx) - c J J ^-V t - S f(y,x) ^_a(dy)ds], 
where Vtf is given by the unique solution of the nonlinear p.d.e. 

Vo/ = /£C +(7x(0,co)). 

Therefore there is a unique limit point and the proof is complete. ■ 

For fixed c > 0, a > and arbitrary atomless measure a on I = [0, 1] with 
total mass a, let £,i(t, dy, dx) be as in Proposition 13 . 2 . H Consider the marginal 
process 

t(t,dx) := / ti(t,dy,dx). (3.2.7) 



We call a super subcritical Feller branching diffusion with individual immi- 
gration (SSFBI) with initial state n(dx) = Jj fj,o(dy,dx) (and parameters a, c). 
The expression for the Laplace functional (|3.2.5H (with / only a function of 
x) shows that this coincides with the so-called immigration process with immi- 
gration rate ac corresponding to the entrance law n t (given by (|6.2.12|l in the 
Appendix) and starting from zero measure at time 0. The existence of a super- 
process with immigration corresponding to an entrance law was first established 
by |LSj (Thm. 1.1). The resulting Laplace transform of £(t) with zero initial 
measure is given by 

E exp (-(£(*),/» =exp(-acj K 0+ (V s f)ds^ , /ed(0,oo), (3.2.8) 

see H3.1.18L 

Corollary 3.2.2 (a) The random measure 



/ x£i(t, dy, dx) = ^ b k (t)8 y 

J(0.oo) 



'(0,oo) 

is a purely atomic finite random measure on I in which the atoms bk{t)5 yk 
correspond to the aggregated mass at time t coming from immigrants of family 
type yk&I- 

(b) For the corresponding family of stationary processes {fjj(t)} t ^WL, the random 
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measures rjj (0) converge to the equilibrium for the process with generator (5 
given by J^S.2.2]) . The equilibrium random measure, £ a has Laplace functional 



Eexp (-(£», /))= exp -ac / K 0+ (V s f)ds , /gCi(0,oo). (3.2.9) 



Proof, (a) The random measure /3(t, <iy) on 7 has independent increments and 
no fixed atoms and is therefore purely atomic (see [K], Chapt. 7). 
(b) Given to < 0, r?|(0) can be decomposed into two parts - one coming from 
the initial value at t and one from the immigration in the interval (i , 0). From 
Proposition 13.2.11 it follows that the immigration parts converge (in the sense 
of weak convergence of probability measures on the space of cadlag functions 
D([to, 0], M 1 (J x (0, oo))) to the diffusion limit with immigration, that is the 
process with generator 0. Next note that for each s > the contribution to 
the aggregated measure at time from the state at time to is stochastically de- 
creasing to zero due to the subcriticality and the contribution from immigration 
on (to,0) is stochastically increasing as to I — oo. Moreover, using the moment 
bounds from Lemma T3. 2. 41 it follows that the family of random measures fjf(0) 
is tight. Therefore we have convergence to £/(0), the equilibrium state for the 
process with generator 25. The representation for the Laplace functional of £ a 
follows by letting t — > oo in (|3.2.8|l . ■ 

Remark 3.2.3 1. Corollaru VS .2.2\ (a) implies that for all 5 > 0, asymptotically 
as e — > only a finite number of immigrant families contribute all but 5 of the 
mass. Each atom corresponds to an excursion from zero for the SSFBI process 
and consists of descendants of only one immigrant family. In fact, we will see 
that asymptotically at the particle level each immigrant family corresponds to 
the descendants of one immigrating particle. 

2. Note that the assumption ?7f(to) M £ M (I x (0, oo)) in Proposition 
\3.2.1\ puts constraints not only on the aggregated mass but also on the family 
structure of the population. To understand what happens if this condition is 
not satisfied consider rjj(t a ) = J2 a k^x E k with J2k a k x k = ® but inf f. x £ k — > oo 
as e — > 0. In this case at times t > to, ??f(t) — ► due to ultimate extinction 
of the critical family level branching. This observation is used below to prove 
by contradiction that the equilibrium populations in Bj, N ^ are asymptotically 
composed of families of size 0(N e / 2 ). 

3. Similarly, if the immigration mechanism is such that it feeds a few large 
families rather than giving small new families a chance, then in the time sta- 
tionary process the family branching makes everything extinct as e — ► 0. 

3.2.2 Moments 

The following lemma was used in the proof of Corollary 13 . 2 . 21 and will also be 
needed below. For the ease of notation we put to = 0, otherwise we would have 
to replace t by t — to- 
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Lemma 3.2.4 For t > 0, let rj e (t,dx) := Jj r)j(t, dy, dx), where rjj is as in 

nnn\) . Let 

m 3 . k {t) = E[(r] e (t),x j ) k ] , j,k£N, 

M(t)=E((v E (t),x)( V %t),x 2 )), 

and o e (l) denote a term that is uniformly bounded in e and converges to as 
e — > 0, and 5(t) denote a term that is uniformly bounded in t > 0, and 5(t) 
converges exponentially fast to as t —> oo, and o(i)| < const -t for small t > 0. 
TTien 

m 1A (t) = a(l - e- ct ) + e- c *mi,i(0), 

TO2 = ±[i - 2e~ ct + e- 2ct ] + m 2 1 (0)e~ 2c * + I^lilM ( e -c* _ e -2c t) 
2c c 

+^(l-e-**), 



(b) 



m 3 ,i(t) = ^ + o e (l) + o(t), 
m4,i(t) = ^ + o E (l) + o(t), 



mi, a (t) = m lj2 (0)e- 2c * 

+ {1 - 4e- ct + 2cie~ c ' + 3 e - 2ct j 
4c 2 L J 

+a 2 {l-2e- ct + e- 2ct } 

+ mi i i(0) _ cte _ ct + 2ac 2 e -ct _ e -2ct _ 2ac 2 e -2 Ci | 

+ m M (0) e 2 o(l -e"** )} 

c 2 

+ — (3a - 2e~ ct + 2ate" 2ct + 4 e - ct (-a + m 1 AO)) + ae~ 2ct - 4mi i(0)e~ 2c *), 
4c 

m 2 , 2 (t) = ^A+ e(^)+5(t), 
2 

M(t)= 4c^ + ^ + ° e(1) + 6W ' 
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Proof. The proof is obtained by applying the martingale problem with the 
generator given by (|3.2.3|) and (|3.2.4H to functions of the form F(fi) — f((fi, <j>)) 
or F(n) = /((/it, </>i), (fi, 4>2)) to derive the following moment equations: 

= ca ~ cmi t i(t), 



dt 

dm 2 ,i{t) 



= TOi,i(i) — 2cm 2j \(t) + cae, 



dm 3 ,i(t ) 
dt 



dt 

cae 2 + 3m2 i(t) — 3cm3 i(t) + 01(e), 



dt 

dmi )2 (i) 



dnii i(t) o 

= cae + 67713,1(4) — 4cm4,i(t) + 01(e), 



dt 
dM{t) 



— 7712, i(t) — 2cm\ t2 {t) + (2ca + e)mi,i(t) + e ca, 



dt 

dmi iS (t) 



cam 2 ,iW + m 3 ,i(t) + m ia {t) - 3cM(t) + 0\{e), 
m^i(t) -4:cm 2i2 {t) +o 1 (e) 



dt 

dm 2i2 {t) 



dt 



= 3cami i2 (t) + 3M(t) — Sennit) + 0\(e). 



Note that the coefficients of the o\ (e) terms only contain moments lower in the 
hierarchy and hence are asymptotically negligible. The results were obtained 
by solving this linear system using MAPLE. ■ 

Remark 3.2.5 In the case in which we replace the constant immigration rate 
a by a random function of time a(-) the expression for m\ l2 (t) becomes 

mi, 2 (t) = mi, 2 (0)e- 2ct (3.2.10) 
, n* u (0) 



{ e -ct _ cte ~ct + 2c 2 ae -ct _ e -2ot _ 2ac 2 e -2ct} 
t pt ps 2 

k 1 (t, s)a(s)ds + / / k 2 (t, s 2 , si)a(si)a(s 2 )dsids 2 
Jo Jo 

t 



+o e (l) • / k 3 (t, s)a(s)ds 
Jo 

where ki(t, •), i = 1,2, 3 are bounded non-negative kernels satisfying 

pt pt ps 2 

sup / ki(t, s)ds < 00, i = 1,3, sup / / k 2 (t, s 2l s{)ds\ds 2 < 00 (3.2.11) 

t Jo t Jn Jo 



t Jo t Jo Jo 

and o e (l) ^ as e — > 
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3.2.3 SPDE representation 

Let £ be an SSFBI process starting from zero measure at time as in subsec- 
tion recall that the Laplace transform of £(t) is given by (|3.2.8[) . By an 
argument similar to that of |LS| (Thm 1.2) it follows that there is a unique 
orthogonal martingale measure M(dsdx) on [0, oo) x (0, oo) having quadratic 
variation measure {M)(ds dx) = ds£(s,dx) such that 

«(*),/> -<£(0),/> (3.2.12) 
' Ms),G c f) + acK 0+ (f)]ds+ [ ( f(x)M(dsdx), / e D P (G C ). 

lo JO J(0,oo) 



Proceeding as in the proof of Theorem 1.7. of LS and tracing the arguments 
of |KS| one infers that has absolutely continuous states, that is, £(t,dx) = 
x)dx, and that one can define a time-space white noise Wt(x) on an extension 
of the original probability space such that 



M(dsdx) — y/£(s,x) W s (x)ds dx. 



(3.2.13) 



Moreover, £(t, x) is almost surely jointly continuous in (t, x) 6 [0, oo) x (0, oo) 
and is a solution of the SPDE l|3.2.18[l below. Note however that in contrast to 
|LS| . £,{t,x) does not have a finite limit as a; J, 0. Indeed, putting 



ft />oo 

z t(x) = / p t -s(y,x)M(ds,dy) 
Jo Jo 



(3.2.14) 



where p is the transition density of the c-FBD process (see section 16.21 in the 
Appendix), we obtain as in LS^ ((4.10), (4.11)): 



£(t, x) = Zt (x) + / acK s (x)ds, x > 0. 



(3.2.15) 



From (|6.2.16|l we obtain 



, . , 2ac / —lex 

acK s (x)ds = exp I — - 

o x \ 1 — e c 



(3.2.16) 



Thus, E(/ 1 £(t, x)dx) — C(t,e)\\oge\ where C(t,e) is uniformly bounded 
away from 0, and 



Var / x)dx 



f(s,y)M(ds,dy) 



o jo 

t />oo 



f(s,y) 2 £{s,y)dy 
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where 



y cic ( — ^cy \ 

Pt-s(y,x)dx < - A 1, and E£(s,y) < — exp I — 



where we have used a stochastic Fubini theorem (cf. jlWj Chapt. 3, Lemma 
4.1). 

Therefore 

Var ( [ £(t,x)dx) < C\{t) + C 2 (t)| loge\ 



for some positive constants C\ (t) , C<i (t) . We then obtain for any 6 > 
P 



N Ci(t)+C 2 (t)|loge| 



Z(t,x)dx-C(t)\]oge\ > 5C(t)\loge\ j < ( ^, (/) ^ 2 

(3.2.17) 

and therefore £(t, x)dx converges in probability to oo as e — > 0. Since 

/ £(t,x)dx is monotone in s this convergence must be a.s., and thus £(t,x) 
is a.s. unbounded as x — ■> 0. 

Finally, recalling (|3.1.13|l and l|3.2.13|l . we see that (|3.2.12|) is the integral 
form of the SPDE 

d 



^x) = y/£(t,x)Wt(x) + G* c £(t, x) - caS' (x). (3.2.18) 

Let us comment on the meaning of the three terms on the right hand side of 
(|3.2.18|l viewed as the limiting family size process as N — > oo in a ball Bi N '. 
The first one comes from the family branching, the second one incorporates 
the individual branching and individual emigration at rate c — eg from B\ N ^ 
(recall that G c is the generator of a c-subcritical Feller branching diffusion), 
and the third term describes immigration of small families into i?« from the 
surrounding medium at a large rate. In fact, S' can be viewed as the limit 
as e — > of a large number 1/e of small families of size e; note that (l/e)6 e 
converges to —S' in the sense of Schwartz distributions on the smooth functions 
vanishing at 0. 

Weak uniqueness of the solution of H3.2.18|l follows from Proposition 13.2.11 
however, as in the case of |KS| it is an open question whether strong uniqueness 
holds. 

The total (or aggregated) population size 



C(t) = / x£(t,x)dx (3.2.19) 

7(0, oo) 



solves the equation 



<K{t) = \ / x)dx dW t - c(C(t) - a) dt, (3.2.20) 

which is a one-dimensional projection of equation (|3.2.18|1 . Note that the process 
£ is not Markov. 
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3.2.4 Equilibrium canonical moments 

As t — > oo, the SSFBI process £(i) with parameters c and a converges in distri- 
bution to the infinitely divisible equilibrium random measure £ a , cf. Corollary 
l3~2~2l Writing 



( a ~ / x£ a (dx) (3.2.21) 

J(0,oo) 

for the "aggregation" of £ a , and v c for the canonical measure of £ a , we obtain 
from l{3~2~%j) and l|3.1.18|l : 

poo 

Vc{[bxM) ■= c / m([b 1 ,b 2 })K t {dm)dt, < h < b 2 , (3.2.22) 
Jo 

According to Lemma 13.2.41 £ a has first and second moments 

E[C a ] = a, E[(C a ) 2 ] = a (a + . (3.2.23) 

Definition 3.2.6 Let us write v c for the size-biasing of v c ( cf subsection TKT]) . 
Because of the well-known relations (cf. Remark |4. 5.2(1 

E[C°]=/ xv c {dx), E[(C a ) 2 ] -E[C a ] (E[C Q ]+ / xv c {dx)\ 

J(0,oo) \ J(0,oo) / 



we obtain immediately from (|3.2.23|) : 
Remark 3.2.7 a) J, Q ^ xv c (dx) = 1 

Remark 3.2.8 We also note that infinitely many immigrant families contribute 
to £ a . This follows from 

v c (0, oo) = oo. (3.2.24) 
To see this note that from \3.2.2ty and yj.l.lty . 

u c (Q,oo) = c (-logP(£° = 0))dt. 



For S > 0, x > let f s (x) = § A 1. Recalling \3.1.T[ and \3.1.1Q)) note that 
P(£° = 0)= lim lime^W^^/ 5 . 

6>— >oo e^O 

T/ien 6t/ a simple modification of \DH$ (6.10). for any S > and > 1 

lim ^)( £ ) >e - tlim ^M£) >e -, r 

o £ " — o e " Js tK ' 
Therefore by \6.2.1b)) 

I* oo />oo 

^ c (0, oo) > clim / e~*( / K t (x)dx)dt = oo. 
5 1° Jo J5 
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4 The genealogy of jumps in a cascade of sub- 
ordinators 



In this section we will carry out the program outlined in Subsection 12.1.11 to 
obtain a representation for the sequence {Q } in terms of a cascade of subordi- 
nators and then use this representation to obtain a genealogical description of 
the population. 

4.1 Propagation of equilibria 

In the preceding section we encountered the equilibrium distribution for an R + - 
valued process whose dynamics is given by ((3. 2. 20(1 . A simpler situation 
is the one corresponding to one-level branching where we have the equilibrium, 

C Q of 

d((t) = y/c^)dW(t) - c(C(t) - a)dt, (4.1.1) 

recall subsection 12 .11 In this section we will derive structural results which are 
common to both situations. We therefore denote the equilibrium states of l|4. 1.1(1 
and l(3.2.20|) by the same symbol £ a . 

In both situations the dynamics has two parameters a and c, and the equi- 
librium distribution is infinitely divisible with expectation a. Therefore as in 
(|2.1.1|l this equilibrium distribution has a representation as 

C(C) = C(S(a)), (4.1.2) 

where S(t),t > 0, is a subordinator with E5(r) = r, r > 0. 

Let us denote the Levy measure of S by fi, and note that /i has expectation 
1. Note that fi = v c (defined in subsection I3.2.4|) if ( follows the dynamics 
(|3.2.2L)[1 . and \i — j c (given by l|t).2.14ll % l if ( follows the dynamics (|4.1.1|l . 

S(a) has a Levy-Khinchin representation as a Poissonian superposition 

S(a) = d W. (4-1.3) 

i:Ti<a 

where '^2 i S( Ti ,yi) is a Poisson population on M + x (0, oc) with intensity measure 
dTn(dy). Since by l|6.2.14|l . resp. (|3.2.24|l . 

H(0, oo) = oo (4.1-4) 

in the two cases, P(S(a) — 0) = if a > and S(-) has infinitely many jumps 
in any open interval. 

Since C° is a Poisson superposition of immigrant clusters (recall 1(3.2.22(1 
and ((6.2.13(1 1. the representation 1(4.1.3(1 has a genealogical interpretation: the 
summands yi measure the size of those parts of £ a which trace back to one and 
the same immigrant. 

Our aim in this section is to study the hierarchy 1(2.1.11) into which 1(4.1.11) 
and 1(3.2.20(1 are embedded. In both situations, the parameters of the hierarchy 
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are a sequence (ci)i=\ t 2,... of positive numbers. Recall that, for different levels 
I in the hierarchy, the dynamics of Q(t) run at separated time scales, and 
the equilibrium state (e+i at level I + 1 acts as (random) parameter a for the 
dynamics at level t. Instead of one Levy measure fi, we now have a sequence of 
Levy measures (fie) ( which is cither (j Ct ) or (v Cl )). 

In this way we obtain a Markovian dynamics which transports the equilibria 
down the levels: 

given = a , Q is infinitely divisible with canonical measure apu. (4.1.5) 

Since \i has expectation 1, (<^) constitutes a backward martingale. We now turn 
to the following problems: 

a) Find a condition on (jii) which guarantees the existence of an entrance 
law, denoted by for (Q) starting in 9 > "at level oo" and having constant 
expectation 9. 

b) Describe the "branching genealogy" underlying such an entrance law. 
We will answer these questions in the next subsections. Later on, we will give 

a relation with the asymptotics of the genealogy of the equilibrium branching 
populations on as N — > oo. 

4.2 An entrance law from infinity 

Let Sk, k — 1,2,... be independent subordinators with Levy measures fik- We 
denote the second moment of jik by . For j > £ define the random variables 

Si(a) :=S e (S i+1 (...(S j . 1 (a)))), (4.2.1) 

and write ILj (a, .) for the distribution of Se(a), I = 1,2,..., a > 0. 

Proposition 4.2.1 If 

oo 

^m k <oo (4.2.2) 
fe=i 

then, for each 9 > 0, the sequence of processes 

(S'^O), Sj_ 2 (6),..., Si(6),Si(0)) 
converges as j — > oo (component-wise) in probability to a sequence 



which obeys 
and 



Se((! +1 ) = $ a.s. for all £ (4.2.3) 

lim <f = a.s. (4.2.4) 

In particular, the distributions ire = £(C* ) are an entrance law for the backward 
Markov chain with probability transition function 

P(QeA\Q + i=a) = U e (a,A), (4.2.5) 
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and they are its unique entrance law with the property 

tt( ^ 5 e as £ — > oo. (4.2.6) 
Proof. Since Sk{a) is infinitely divisible with canonical measure a[i, we have 

Var(5fe(a)) = am k . 

Hence we obtain 

Var(S£ +2 (a)) = Var[E[5 fc (5 fe+ i(a))|S fc+ i(a)]] + E[Var[5 fc (5 fc+ i(a))|5 fe+ i(a)]] 
= a(m k+ i +m k ). 

In the same way we get for all j > k > £: 

Var{S J k (a)) = a(mj-i + ... + m k ) (4.2.7) 

and 

E(Sj(a) - S^{a)) 2 = a(mj-_i + ... + m k ). 

Thus, due to 1(4.2.2(1 . for fixed £ the sequence (S J e (9))j > i is Cauchy in L 2 . We 
define 

C| = L 2 - lim Sf(0). (4.2.8) 
Since a i— » St(a) is continuous in L , we have 

= lim CiW) = lim = C? a.s., (4.2.9) 

j — >oo j— »oo 

which proves (|4.2.Hfi and, a fortiori, implies 1)4.2.5)1 . From ()4.2.8)1 and 1)4.2.7)1 it 

is clear that 

oo 

VarQ 9 = m k . 

k=e 

Since E£| = 0, this together with 1)4.2.2(1 implies that converges to 9 in 
probability as £ — » oo. Moreover, since because of 1(4. 2. 3|) £| is a backwards 
martingale, this convergence is even a.s., and we have 14.2.411 . 

It remains to show the claimed uniqueness statement. For this let (wk) be 
an entrance law for (II^) obeying 1(4.2.6(1 . and let X k , k = 1,2,... be random 
variables, independent of the subordinators St, with 

From the entrance law property of (jr k ) and the definition of (11^) we have for 
all j > £: 

■K l = C{S i l ,{X j )). (4.2.10) 
On the other hand we have by monotonicity of r i— > Sj(r): 

ns{{x J )-si{9)\<nx J -9\. 
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From this, the claimed identity ng — £(£|) follows by Markov's inequality to- 
gether with jSQI and (|4.2.10() . ■ 

Proof of Proposition l2Ti . II 

Recall that with the notation of Proposition 12 . 1 . ll Tig (a, .) = C(S^\a)) and 
11^ (a,.) = C(S < j, 2 \a)) 1 where and S^p are subordinators with Levy mea- 
sures 7£ and i>n, respectively. From Remark 16.2.31 we have that the second 
moment of 7^ equals l/(2cfc), and Remark 13 . 2 . 71 shows that the second moment 
of Vk equals l/(4c|). The proof of Proposition 12.1.11 is thus immediate from 
Proposition 14.2.11 

4.3 The genealogy of jumps in an iteration of subordina- 
tors 

The composition of subordinators gives rise to a "genealogy" of their jumps. To 
illustrate this, consider the two subordinators Si, S 2 , where 

S 1 (b) = J2(Si(U)-S 1 (t i -)), 
S 2 (a) = E - S ^ T n~))- 

T n < a 

Then 

S 1 (S 2 (a))=J2 E (Si(U) - S^U-)). 

r n <a S 2 (r„-)<t I <S 2 (T„) 

In this way, the jumps of S\ are coagulated into families of jumps stemming 
from one and the same jump of S 2 . 

Iterating this, we obtain from the flow property (|4.2.3|l that 

Qt = E {Si{ n )-Si{n-)) & s. , (4.3.1) 

n£[0,Cf] 

where {t^} is the set of all points in [0, Cj) in which t i— > Sj(r) has a jump. 
The representation (|4.3.1|l induces a partition of [0, Q] which we denote by 
y($j,t. Note that for fixed £ the tyj.t are coalescing (i.e. becoming coarser) as j 
increases. 

The sequence of coalescing partitions tyj.i induces a graph (S^ as follows: 
The set of nodes of <S>£ is the union Uj>^{i} x For n ^ &i we call its first 
component the level of n. For two nodes n\ — n 2 = (j 2 ,I 2 ) of <S( we 

say that ni is an ancestor of n 2 if ji > j 2 and I 2 Q lit an d we say that n\ is the 
•parent of n 2 if ri\ is the ancestor of n 2 with j\ = j 2 + 1. The (directed) edges of 
35^ then are all the parent-child pairs in ®g x &£. Say that two nodes in are 
related if they have a common ancestor. Then by construction of the sequence 
PPj,^) each equivalence class of &g is a tree, i.e. a directed connected graph 
without cycles. Therefore, <&e is a forest, i.e. a union of pairwise disconnected 
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trees. Finally, we label each node of <&i with the length of the subinterval of 
[0, Q] to which it corresponds, thus arriving at the random labelled forest 
which we associate with the random sequence of coalescing partitions 9pj,£, and 
which encodes the genealogy of the jumps of the process (..., Q) constructed 
in Proposition 14. 2 . ll 

Intuitively, viewing [0, q] as a continuum of individuals, this means that 
two individuals a±, a-2 G [0, Q] belong to the same element of tyijj if and only if 
they descend from a common ancestor (or equivalently, from one subordinator 
jump) at some level less or equal than j. Furthermore two individuals a\, a% G 
[0, (%] belong to the same element of the minimal partition ^oo,^ if and only if 
they descend from a common ancestor at any level higher than t. Using the 
independent increments property in 6, l|4.1.4l) and (|4.2.9|l it can be shown that 
there are countably many distinct elements in ^oo^ each corresponding to an 
infinite tree. Therefore we have a decomposition of the equilibrium population 
into a countable set of subpopulations each consisting of individuals having a 
common ancestor. 



4.4 The genealogy in the hierarchichal mean field limit 

With the special choice of (Si) described at the beginning of this section, we 
have all reasons to conjecture that the random labelled forests $e defined in the 
previous subsection describe the genealogy of the (one- or two-level) branching 
population in equilibrium as N — > oo. 

To make this more precise, consider a fixed sequence B^ N \ £ — 1,2, ... of 
nested balls in fijv ; and let Vf be that part of the equilibrium population 
which lives in Bg N ^ . (Here and below we suppress the notation of 9 > which 
we keep fixed.) Fix N and £ for the moment. For two individuals I\, 1% in V^ N ^ 
and j > £ we say that 

h ~i h 

if I\ and I2 have a common ancestor in Bj N ^ . This induces a partition on v\ 
which we denote by ffSji- The sequence 9PjP, j = £,£ + 1, ■■■ is coalescing, 

and we can associate with it a labelled forest in the same way as we 

associated $e with tyj.e in the previous subsection, the only difference being 
that now we label the nodes of by the cardinalities of their corresponding 
sub-populations of V t , divided by N l . 

Our main result (Theorem 12.4. lfl suggests to conjecture that, in a suitable 
topology, 

#° -> Si as N -> 00. 

4.5 Size-biasing iterated subordinators 

With a view towards the genealogy of a sampled individual (see subsection 14.6(1 
we will now prove a representation of the size-biasing of C(Sj(a)), where Sj is 
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the composition of subordinators defined in 14.2.1|) . To this purpose we first 
consider a single subordinator S evaluated at a random argument. 

Proposition 4.5.1 Let S(t),t > 0, be a subordinator with Levy measure de- 
noted by jj,, and let A be an R + -valued random variable independent of S and 
with finite expectation. Then the size-biasing of C(S(A)) arises as the distribu- 
tion of S(A)+Y, where C(A) is the size-biasing of C(A), C(Y) is the size-biasing 
of fj,, and A and Y are independent. 

Proof. We write II rM for the distribution of a Poisson point configuration on 
K + with intensity measure r\i. Then 

£(S(r))=£ T «*,id R+ )), 

where 

£ r (*)=n r(1 . 

Writing a for the distribution of A, we thus have 

C(S(A)) id R+ )), (4.5.1) 

where 

= J n TM (.)a(dr). 

Our task is to compute the size-biasing of £ CT ((\]/, id^,)) with s(ip) = (tp , idn + ) , 
cf. Definition I6.1.1I To this end let us first compute the size-biasing of At ('I') 
with (t/j, idi£ + }, and then project. It follows from Corollarv lG. 1 ,4l in the Appendix 
that the size-biasing of £ CT (4') with id]R + ) is £($ + 5y ), where 

£(*) = J U T ^.)a(dr), 

a is the size-biasing of a, C(Y) is the size-biasing of fi, and $ and Y are inde- 
pendent. Consequently, the size-biasing of £ CT (( 1 J r , idfi + )) with (i/>,idR + ) is 

£((<Z> + Sy,id R+ ))=£(S(A) + Y)), 

where C(A) = a, and S, A and Y are independent. Together with l|4.5.1|) this 
proves the claim. ■ 

Remark 4.5.2 For deterministic A, Proposition \4-5.T\ renders the well-known 
fact that the size-biasing of an infinitely divisible distribution 7r on R + is the 
convolution of n with the size-biasing of the canonical measure ofir. 

Corollary 4.5.3 Let Sj(a) be the iteration of subordinators defined in 

where the Sk are independent subordinators with Levy measures fik- Then the 

size-biasing of C(Sj(a)) arises as the distribution of Si (a) defined by 

Sj(a) = Si(a) + Sr\Y 3 -i) + Sf 2 (Y^ 2 ) + ... + + Y e , (4.5.2) 

where C(Yk) is the size-biasing of [i^, S\ is distributed as S\ and all random 
variables occurring on the r.h.s. of {^.5.^ are independent. 
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Proof. For^ = Proposition ^. 5. ll shows that the size-biasing of C(Sj 1 (a)) = 
C(Sj-i(a)) arises as the distribution of 

Sj-x (a) +y J _ 1> 

where both summands are independent. One more application of Proposition 
I4.5.1l thus gives that the size-biasing oi £(Sj_ 2 (a)) — £(5 7 -_2(S'j_ 1 (a)) arises as 
the distribution of 

V2(Si-i(a)+Y--i)+%-2 
which due to the independence of Sj-i(a) and equals in distribution to 

Sj-iiSj-ila)) + S'^iY^x) + Yj- 2 , 

Sj_ 2 being an independent copy of Sj-2- 

Iterating the argument we arrive at our assertion. ■ 

Remark 4.5.4 As before, let us denote the second moment of jik (or equiva- 
lently the first moment of Y^ ) by rrik ■ From [4.5.^) it follows that 

i-i 

ESj{a) =a + ^2m k (4.5.3) 

k=e 

Hence the summability of the is a sufficient condition for tightness of the 
SI (a). 

We now turn to the entrance law constructed in Proposition I4.2.T1 

Corollary 4.5.5 From {4-2.8\j and Corollary \4-5.S\ we obtain that the size- 
biasing of £(Cf) arises as the distribution of 

Cf = Y e + S e +\Y e+1 ) + S f ( + 2 (Y e+2 ) + .... + & (4.5.4) 

where the random variables Yk and are as in Corollary \4-5.2\ and all random 
variables occurring on the right hand side of {4-5. 4\l ar ^ independent. 

We can go one step further and study the genealogical relationships under- 
lying the representation 14.5.4|) . 

To this purpose, let us study the branching dynamics on the "populations 
of jumps" induced by the composition of the subordinators Sk, resuming the 
reasoning of subsection l4.3l For each k we consider a branching dynamics which 
takes a counting measure <fik on R + into a random counting measure $fc_i in 
the following way: If <pk — Sie/ fc then 



ieik 



where Vf^ is a Poisson counting measure on R + with intensity measure yi^—i, 
and the are independent. 
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Now fix two levels j, £ G N with j > £. Starting with <j)j = 5 a , and iterating 
the branching dynamics from level j down to level I, we construct a random 
path of counting measures on R + which we denote by 

(S a ,^ j _ 1 (a),...,^i(a))=:Hj{a). 

By keeping track which atom in < &^_ 1 (a) stems from which atom in ${(a), we 
can enrich the history Hj (a) to a tree TJ? (a), each of whose nodes is marked by 
a non-negative real number. For example, if $^(a) = Ylj^j S yi , j > k > t, then 

the set of nodes of Tl (a) at level k corresponds to the index set Ik , and yi is 
the mark (or "size") of the node with index i. We write 

fi(T/(a)) = ($Ka),id R+ ) 

for the total size of the tree (a) at level I, and note that 

*(T/(a)) = rf Aj(a), 

where Sj (a) is defined in H4.2.1|) . Proceeding in a similar way as in the proof 
of Corollary 14.5.31 we obtain a "spinal decomposition" of the size-biased tree, 
which we state here without proof. 

Proposition 4.5.6 The size-biasing of £(Tj!(a)) with s(Tj!(a)) arises as the 
distribution of the superposition of Tj? (a) and 

(s^,. . . , Sy k , . . . , k = j-i,j-2,...,e, 

where £(Yjt) is the size-biasing of fj,k, Tjl (a), Yj—i, . . . , Yt are independent, and 
given Yj-i, ■ • ■ , Ye, the $"(1^), k > r, are independent. 

Under the assumption (|4.2.2[) of summability of the second moments of Hk, 
HN, the size-biasing of the random labelled forest 3^ = $£ (constructed in 
subsection 14. 3|) with respect to its " size" s(g^) = £f arises as the independent 
superposition of $f and j^°°' can ; where j^°°' can i s constructed as follows: 

First build a "spine" (. . . , Y^i, Xj-i) Yf), and given the spine, superimpose 
independently the trees T^(Yfc), k > I. 

4.6 The genealogy of relatives of a sampled individual 

In subsection 14.31 we fixed a ball (or i?-block) Bg from the beginning. Now 
we take a different viewpoint and think of an individual sampled from the equi- 
librium population within a union of many ^-blocks in Qn from the beginning. 
Denote the chosen individual by /, and the i?-block by Bp. Recall that for large 
N the total number of individuals in an ^-block is approximately distributed like 
N Q (see Theorem 12.4. 1\ . we see that the number of individuals in the chosen 
block is approximately distributed like N l Q, where £{Q) is the size-biasing of 
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The block Bg sits in a nested sequence of blocks of levels £+1,1+2, . . . which 
we denote by Bg + i, Bg + 2, . ■ . Consider the population Vg of all those individuals 
in £>£ which have an ancestral family in common with the individual /. The 
population Vt can be decomposed in a natural way according to its immigration 
history into the Bj, j > £. 

For j > £, denote by Vl the subpopulation of all those individuals in Vt 
that have some common ancestor with I who lived in Bj but none who lived in 
Bj—\. In other words, V\ consists of all those individuals J which obey J ~j I 
but not J /. 

In this way we obtain a decomposition of Vi according to the hierarchical 
distance of the (geographically) closest ancestors common with the chosen in- 
dividual /: for j > £, the subpopulation V\ consists of those individuals in Bp 
whose geographically closest common ancestor with / has hierarchical distance 
j from /. 

The size (i.e. the total number of individuals) of V\ is approximately dis- 
tributed as 

N l S\-\Y 3 - X ) = N l St(S t -x ■ ■ ■ (Si_2&-i))) 
and thus has approximate expectation N -r-k — ■ Hence the summability con- 

dition (|4.2.2|l (which corresponds to the condition for transience resp. strong 
transience of the hierarchical random walk) amounts precisely to an expected 
finite number of relatives of the chosen individual in the block Bg. 

5 The hierarchical mean field limit of two-level 
branching systems in equilibrium 

In this section we investigate the two level branching equilibrium ty( N ' S \0) 
described in Proposition 12 . 3 . 51 and its limiting behavior as N — > oo. We assume 
that the underlying random walk is a (2, (q), A^)-random walk on f2jv and (q) 
satisfies the strong transience conditions (|2.3.9|l and (|2. 3.10(1 . Recall that (B^) 
denotes a sequence of nested blocks in fijy. 

We will see in Lemma 15.2.11 that in equilibrium, asymptotically as N —> oo, 
$( JV ' e )(0) consists of the order of N e / 2 families in B^ N \ a typical such family 
having a random multiple of N e / 2 individuals. 

5.1 A spatial ergodic theorem 

In this subsection we first collect some basic facts about the two level branch- 
ing systems Sfr( N \t) on which were introduced in section 2.2.6 and their 
equilibria. The main result will be a spatial ergodic theorem for the aggregated 
equilibria; this will also be an ingredient in the proof of Theorem l2.4.1l 

Let M c (f2jv) \ {o} denote the space of non-zero counting measures on 
such that finite sets have finite measure, and let M 1 {M (f2jv) \ {o}) be the set 
of measures v on Mc(fijv) \ {o} such that J M (q \\r ^^(B)v{d^i) < oo if B 
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is a finite set. Then let M^(M c (Qn) \ {o}) denote the subspace of counting 
measures in M 1 (M c (n^) \ {o}). The set M c (£In) \ {o} carries a complete 
separable metric p which generates the restriction of the vague topology on 
M c (fl^) to M c (17jv) \{°} and makes a subset of Af c (17jv) bounded if and only if 
it is contained in {[i\n(B) > for some finite B C £In}, see MKM Proposition 
3.3.2. The set M 1 (M c (fijv) \ {o}) is equipped with the topology generated 
by is i ► J M (n„)\{ } M-^OK^m); where B is a finite subset of fijv, and v i— > 
J F(/j,)i'(dfi), where F : M c (£1tv) \ {o} ^ K is continuous, bounded and with 
p-bounded support. 

We can represent ty( N >(t) as a cadlag Markov process with state space 
M c 1 (M c (fi A r) \ {o}). On M^(M C (Q N ) \ {o}), we consider the class of func- 
tions D(<S nN ) = {H{v) = h (J F(p)v(dn))} with = /((/x,^)) where <p 
has finite support and h, f are continuous functions on R with bounded second 
derivatives. We also define 



Sp(x) de 



F{fi + eS x )\ £=0 , 



t) 2 



5fx(x)8fj,(y) deide 2 



F(p + e 1 S x +e 2 5 y ) 



£l=£2=0' 



Then *W is the unique solution to the martingale problem given by the 
generator 



<8 Qn H(v) 



(5.1.1) 



J2 J v{dn)h!{ J F(ji)v{dn)) 
5F(fx) 



xen 



E 

\k=l 



Cfe+1 

N k / 2 



(fj, k (x) - (t(x)) 



S 2 F(n) 
5/j,(x)5fi(x) 



where 



h"( / F(»Mdri)F 2 (»Mdn) 



1 



= -^k £ /*(»)• 



yGQ N :d N (y,x)<k 

Proposition ^ .3 . 5l establishes the existence of a non-trivial equilibrium ty( N ' e '> 
that is spatially homogeneous (that is, with law invariant under translations in 
Hn) and mean 0. The equilibrium random field {i/j x N (t)} x ^n N is defined by 



(5.1.2) 



This gives the total number of individuals at site x irrespective of their family 
memberships. 

Recall that {i]^ N ' e \t,dx)} defined in 12.4.11f) describes the equilibrium nor- 
malized family size process in the block Bg N ^ with mean number 6 individuals 
per site. 
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Also recall that the normalized equilibrium population mass in B# is given 

hv rsee Exmrcxra 

cf» e) W = f X r ) f' 6 \t,d x ) = l l ^ fi \N^t) 

= W E f K*)* {N > e) (N" 2 t,dn). (5.1.3) 

Lemma 5.1.1 Lei Gjv denote the Green operator of the (2, (cj), N) -random 
walk with (cj) satisfying condition ]2.3.1(]\) for strong transience, and write, for 

(yi, Va) = E ^1(^)^2(2;). 

TTiera 

{<Pi,G N <p 2 ) < const 2^ ^1^1/2 i 5 - 1 - 4 ) 

00 1 

(<Pi,G%(p2) < const ^i(x)(^ 2 (t/) ^ — (5.1.5) 

where the constants do not depend on N . Hence 

(<pi,G N ip2) < const(l,(^i)(l,</j 2 ), (5.1.6) 
((pi,G 2 N ip 2 ) <const(l,<^ 1 ) (1,(^ 2 ). (5.1.7) 

Proof. The transition probability pt(x,y) of the (2, (cj), iV)-random walk 
is given by 



j=\x-y\ + l 



where the hj N ^ are positive numbers (depending on N) such that 

h {N) > const C3 ~ 1 ? > 1 

ft 3 ^ const nU _ 1)/3 , 

(see |IX;W2p . Hence 

G N ip{x) = y2p t (x,y)ip(y)dt 
Jo y 

OO 

<(N-i)J2v(y) E 
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<comt(N-l)J2<P(v) E -J^— 
y i=|x-i/|+i 3 

oo 1 

< const >(z>(w) > — — - — 

y j=\x-y\ J 

< const E jv*/2 

y j=\x-y\ 



< const ^(y)—- 



and then l|5.1.4|) follows. 
Similarly (see |l)(4W2p 



G%tp(x) = 



J t {YlPt{x 1 y)^{y) S jdt 



<(N-l)^(v) £ 



1 



=t|+i^ J (^) a 



OO A T 1* _ 1 



<const(JV-l)5>(y) £ ^ 

y 3=|x-y|+l 3 1 

oo 1 

< const 2jv(y) 2J ^2 

and (|5.1.5|l follows. ■ 

Remark 5.1.2 In 'DGW21 the exponential waiting time of the hierarchical ran- 
dom walk has parameter 1, whereas here this parameter is J"^ j , see subsec- 
tion 12.2.H In the present case, YleQt — Y^jLo =: Ln- This amounts 

to a time change t — > L^yi m i/ie calculations above, which produces L^h]- in 
place ofhj and does not change the results of Lemma \5.1.1\ ( assuming cq > 0). 
The following is the analogue of [DGW1] (2.3.3). 

Proposition 5.1.3 Under the conditions of Lemma \5.1.1\ the first and second 
moments of ^ N '°\o) are given by 

and 

E (Ci N,e \0)"j = 2 + 0((<PN,t,<PN,l) + (<PN,t, G N ip Ny i) + ^(ip N j,G 2 N ip N j)) 

(5.1.9) 

where tpN,e( x ) = ^^- B ( N A X ) an d Gn is the Green operator of the (strongly 
transient) (2, (eg), N) -random walk. 
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Corollary 5.1.4 If the (2, (cj), N)-random walk satisfies conditions \2.2.2\) and 

\2.3.10\) . then for each i 

(a) 

snvE{{cf' e \0)) 2 ) < oo. (5.1.10) 

N 

(b) 

Vax(C { t N ' e) (0)) ^0 ast^oo (5.1.11) 

uniformly in N . 



Proof, (a) follows immediately from Proposition 15.1.31 and Lemma f5. 1.11 
(b) From Proposition 15.1.31 

Var (Cf ^(0)) = 0(((p N j, ip N ./) + ((p N ,e, G N <PN,i) + -7{<PN,e, G 2 N <j> N ,i)) (5.1.12) 

We will show that each of the three terms on the r.h.s. of (|5.1.12|) converges to 
as I — > oo uniformly in N. First {<fN,e, <PN,i) < -p ^ as ^ - > oo uniformly 
in TV. 

We have from H5.1.4|) and using the ultrametric property of • |, 

(<PN,Z,G N <p N! t) <COnst-^ N \x-y\/2 

\x\,\y\<i 



N 2£ \ ^ N\y\/ 2 

x^y, \x\<\y\<t 



T2l I ^ Z-j N k/2 J 



< const NU 

1 / , ]\f( e + 1 ) 3 / 2 - iv 3 / 2 
= const— - N 



N 21 \ N 3 / 2 - 1 

< const ^ + jlj^j ^Oas^oo 

uniformly in N, and from (|5.1.5|l . again using the ultrametric property, 

1 °° 1 

(<PN,l,G%fN,t) < const ^ Y E ^2 

\x\,\y\<ej = \x-y\ °i 



< const W^g_L + ^ g_L 

\ 3=0 3 x^y,\x\<\y\<tj = \y\ 3 ) 
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/ oo jAt \ 



< const 



const 



_1_ ^ 1 1 ^1_ 



(5.1.13) 



The first term goes to zero as £ — * oo, the second term goes to zero as £ — > oo 
by dominated convergence, and clearly the last term goes to zero as t — > oo. 
Therefore 



(<Pn,£, G 2 N (fN,i) — ► as f — > oo uniformly in iV. 



(5.1.14) 



Remark 5.1.5 (a) 



OO -j 

jim o Var(cf^(0))=^E^- 



c ■ 



(b) 



(1 1 °° 1 \ 



Corollary 5.1.6 Consider the "exterior function" 

oo 



fc=l 



Then 



E 



E ^V A c*t(x) f (0,0^) 



< OO 



uniformly in N . 

Proof. As in Proposition 15.1.31 (from |DGWlj . comment 2.3.5) 



= 2 {1, <PN,£,cxt) 2 + ^((ViVAext) <PN,£,cxt) 
+ (fN,e,cxt,GNfN,£,cxt) + - ( t PN,i,e^t)GN ( PN,e,e^t))- 



(5.1.15) 



(5.1.16) 



(5.1.17) 



(5.1.18) 



(5.1.19) 
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We will prove that each term on the r.h.s. of (|5.1.19|1 is bounded uniformly in 
TV. 



(1, <PN,e,ex.t) = / „ 

fc=l 
oo 



fc=i 





-l 


Ct+k- 


1 


]\[f.+2k- 


-l 


ce+k- 


i 



fc=l fe=l 



< const ^ ' -r /. ■ ■ i 
fe=i 



< oo uniformly in TV. 

Then, from (|5.1.6[) the terms (vJjv,*,ext, Gjf<pN,t,ext) and (ipNj,cxt,G 2 N ipN,e,cxt) 
are bounded uniformly in TV and since <fN,i,ext has a uniform bound in TV as 
well, using 12.2.2(1 we see that (<p n, e, cxt, <PN,e, cxt) is also bounded uniformly in 
TV (with TV > supQ + i/c£). ■ 

Theorem 5.1.7 (Spatial ergodic theorem) The pointwise ergodic theorem 
holds on Qjvj that is, 

lim 6 Nfi \o) = 8 a.s. (5.1.20) 

Proof. First note that the equilibrium random field {i/'i'^Jigsj,, defined by 
(|5.1.2|l is invariant under the action of the group n N , and E(ipi N ' e ) = 6. More- 
over by l|5.1.11|) the spatial averages Q N (0) = frEjgBC) ">Px N ' e \o) satisfy 

lim Var(C (JV ' e) (0)) =0, 

j—KX 

so that the convergence in probability follows. To complete the proof note that 
Q,n is an amenable group and the collection of balls {Bg }^gN is a tempered 



F0lner sequence. The a.s. pointwise convergence then follows by [Linj ( Theorem 
1.2). ■ 



5.2 Reduction to two successive scales 

In the section we show that the analysis of the multiscale behavior can be 
reduced to the case of two successive scales. 



5.2.1 The equilibrium family size process 

Let {(i N ' 8 \t), t 6 M}i=...,2,i denote the collection of equilibrium block average 
processes in the nested sequence of blocks B { p (see |£E39). 
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Lemma 5.2.1 The fixed time marginal distributions {Cg N ' S \ty} a t time t = 
satisfy: 

(a) For each t G N the family {(g N ' e \o)}NeN is tight and every limit point has 
expected value 6. 

(b ) The family-size constraint 



E 



xr)f N (0, dx) 



K 



< 



const 



(5.2.1) 



is satisfied for a suitable constant not depending on K > and N > 2. 



Proof, (a) follows immediately from Corollary 15. 1.41 
(b) Let yf{x) = ^rl Bm (x). Then 



K 

< — E 
" K 



xj]f' e (0, dx) 



< — E 
" K 



x 2 4 N ' 6) (0,dx) 



x 2 4 N > e) (0,dx) 



1 1 
1 



E 



-N e ^ 2 E 



Hf)) 2 fW)(0,^) 
(M^f)) 2 * W) (0,d M ) 



A' 



/ i 

< — ^— const 



< 



K 
const 
~K~' 



1 



(by (2-3.1)) 
(by the proof of Corollarv l5.1.4|) 



Remark 5.2.2 Part (b) of the Lemma \5.2.1\ implies that asymptotically the 
restriction of^ N,s \0) to the ball B^ N ^ consists of families of size 0(N e / 2 ) or 
smaller. 



5.2.2 Distant immigrants 

In order to establish that there is a unique limit law and to identify it we now 
return to the dynamical picture. In the next lemma we show that the expected 
contribution to the equilibrium population in a ball Bi coming from immi- 
grants who immigrate directly from outside the ball B^\ and the descendants 
of the population in the ball B^ in the distant past are both negligible. 
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Lemma 5.2.3 (a) Let {^^^(0)} denote the contribution to the equilibrium 
block average in the ball B^ at time coming from individuals immigrating 



directly from (B^l) c . Then 



(C»0) 



(b) The expected mass to enter B^_\ from (B^\) c in a time interval of length 
N^isO^). 

(c) The expected contribution to the population in the ball B^' at time from 
the descendants of individuals alive at time to is of order 0(e _Cf ' to ') as to — > 
—oo. 

Proof, (a) To verify this, we first note that the total number of individuals to 
immigrate from B^S w ^ n k >2to B^ in the time interval [— N e / 2 t, —N e / 2 (t+ 
dt)) is of the order 

&{-t)N^) x (j^B^) >< M >< **' /a 

where the first factor is the number of particles in B^ k , the second factor is the 

rate of migration of each of these particles to a point chosen randomly in B^) k , 

the third factor is the probability this point falls in the tagged ball B[ N ^ and 
the last factor is the length of the time period. Recalling that the mass process 



in B\ ' is subcritical with parameter , the expected total mass at time 



coming from immigration from outside B { /1\ in the time interval [-N e / 2 T, 0), 



T > 0, is of the order 



k>2 

< const x 



—T 



N 1 / 2 ' 



uniformly in T and £, where we have used the assumption ()2.2.2p . Therefore the 

expected mass (normalized by w) at time in I?j that immigrated during 

the time interval [N e / 2 to, 0) directly from outside b\^\ is 0{N~^) as N — > oo. 
(b) and (c) follow from a first moment calculation. 



5.3 Diffusion limit of the family size process in two spatial 
scales 

In this section we consider the asymptotic (as N — > oo) time development of 
the population occupying B^_\ and in particular the subpopulation obtained by 
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considering individuals that occupy a tagged ball B^' in the natural time scale 



for the population in B^ N ' and assuming that the initial family size processes 



in B^_\ satisfy the family size constraint of Lemma T5. 2. II 

As a result of the previous lemma, asymptotically as N — » oo and to — » — oo 
the equilibrium population in s| consists of the descendants of immigrants 

coming from B^_\ during the time interval (—N~ e / 2 to, 0] (recall Remark 12. 3 .4( 1 . 
Another key point is that, as we verify below, in this time interval the family 
size process in B^J is asymptotically constant. Moreover from Lemma !?). 2. H all 
but an arbitrarily small proportion of the population in B^_\ is structured into 
0(N^ +1 ^ 2 ) families containing 0{N^ +1 ^ 2 ) individuals. 

We next show that in the ./V^/ 2 -time scale the total population structure in 
the ball B^!\ is essentially constant. 

Lemma 5.3.1 Let{^£+f\s)} be the equilibrium normalized process (see \5. ) 
in its natural time scale N^ +1 " 2 . 
For t = t (N) < such that ^ < c, 

P ( sup \C%f\N-VH)-c£f\N-VH )\>K) 

/ const \ , . 
< f-p-)c + o(l), (5.3.1) 

where const does not depend on t and N, and o(l) converges to as N — > oo. 

Proof. Recall that $W(t) is characterized as the unique solution of the mar- 
tingale problem with generator (|5.1.1|l . Applying the generator to the function 

F ( v ) = I (M) NT+r^- n (.N))v(d^) : it follows that {Mt +1 (t)} t() <t<o, defined by 
D i+i — — 

M i+1 (t)= (£f\t)-C^f\t ) 

E^(ciK fe (^)-cffi e) ( S ))^ 
fc=i 

is a martingale. 
Then 

P ( sup \(£f\N-^t) - $f\N-Wto)\ > k) 
\t <t<o J 

± p ( /I I E (<ffi!*(*- fc/a <) - c^M I* > § )) 

V 37172 fc=i / 
+ /'| sup |M m (t)| > * | (r,.3.L>) 
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Recall from Corollary 15. 1.41 that E fcj^'^CO) is bounded uniformly in N 
and t. Let 

oo 

a*) -iE^(cffi 9 i(^ fe/2 «)-ca e) ( S ))i 

fc=l 

and note that by Corollary 15. 1.61 E(gjy(s)) is uniformly bounded in s and N. 
Then 



1 ^(^>y) 

C P7? / 



const Jp_ 
S K 2 ( N^ 2> 

const, to 2 



Afl/2 



to 

- 1 / - 



g N {s)ds 



< 



< 



^2 ^1/2' 

const / io 



*o 

■ i i 



E(g N (s) 2 )ds 



K 2 XN 1 / 2 



Next we note that by Lemma l3.2.4f b) and Remark |3.2.5I fwith s = -k) we 
?ct for t < t < 

OO 

™(rAN,e),.s AN, 9), ^ 2 ^(iV,e)/ s _ V"^ Ce+k f\-nf\ 



k=l 



(m e 



2 -2c e+1 (t-t ) 



m 2 ) 



1 



-c e+1 (t-t ) + q_ to j c 2 ee -«a+i(*-*o)j 



H — 2 — / fci(s,t)a(s)ds + / / fc 2 (^, S2, si)a(si)a(s 2 )c?si(is2 

C ^+l "'to J t J t 



+o(l) • / k 3 (t, s)s(s)ds 



la 



where •) i = 1, 3 and &2(t, •, •) are bounded non-negative kernels that satisfy 
(|3.2.11fl and o(l) converges to as N — > oo. Using Corollary 15 . 1 . 41 we conclude 
that 



We now apply the L 2 -martingale inequality. 



^ C0nst 'Vv^ + Ol( iv :) - (5 - 3 ' 3) 



P ( sup \M t +i{t)\ > |- M e+1 {t /N^ 2 ) = 
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<^e(|m^ 



(0)| 2 |M £+ i(t o /iV 1/2 )=0) 



< 



const | to I 
K 2 N 1 / 2 



Let M 1 [[0, oo ) x (0, oo)]\{0} denote the set of non-zero Radon measures 
on [0, oo) x (0, oo) such that J 0+ J yn>{dx,dy) + J Q yfj,({0}, dy) < oo and 

Jo+ /o°° x Kdx, dy) < oo. 
Proposition 5.3.2 Let I be fixed and 

\ N e/2 ^ ax ' N (£+i)/2 ta «-r 

Assume that 

yW(t ;dx,dy) N ^° ti°(dy)6 (dx). 

Then {N 1 ^ 2 y( N ' > (t;dx,dy)l x> o}t <t<o converges weakly as N — > oo to the 
M 1 (0, oo) -valued diffusion {v(t)}t <t<o with generator (St defined as follows. 
Let F(ne) = f({ne, <p)) with ip G C 2 ((0, oo)) with \ip(x)\ < const\x A 1|. Then 

e t F(fjL e ) =/ / ((w,v»(w,GgVi> + 5/"«W»W»<W,^> (5-3.4) 
+f{{fJie,<pt))ct J ^■(x,y)\ x=0 yn°(dy). 

(2) 

where G CI denotes the application of the operator G ce to the x variable. That 
is, rj(t) is a two-level branching process with constant multitype immigration 
source with immigration from zero of type y at rate given by y/j, (dy) and total 
immigration rate of J y[i°(dy) . 

Proof. The proof is a refinement of the proof of Proposition 13.2.11 We 
begin by noting that Lemma 15.2.31 implies that asymptotically as N — > oo the 
contribution of immigrants into b\^\ from B^f^, k > 2 is negligible. Moreover 
the contribution of immigrants into B^l from B^f^ in the time scale N l l 2 is 
also asymptotically negligible. Therefore in the time scale N l l 2 we can restrict 
attention to the population in B^\ . More precisely, asymptotically as N — > oo 
the population { { e N ' e) (*) = // xN x / 2 y^ (t, dx, dy) lx>o consists of descendants 
of immigrants entering from b\^\ in the time interval [N e / 2 to(N), 0] provided 
that we take to(N) — » — oo. 

As in the proof of Proposition ^ . 2 . ll a standard argument yields the tightness 
of the processes {(e N ' 6 \t)}t„<t<o and {Ce+f\ N ~ 1/2 t)}t <t<a- The tightness 
in C^O^M ^Q.o o))) of N^ 2 y^ N \-)l x>0 is also obtained as in the proof 
of Proposition 13.2.11 One difference is the presence of additional terms in the 
expressions for the moments. However these expressions tend to zero as N — ► oo. 
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A first moment calculation shows that the expected contribution of individu- 
als who leave and then re- immigrate is 0(-^) and therefore asymptotically 

negligible. This means that we can treat the population in Bg N ^ in the time 
interval [to,0] as a two level branching system critical at the family level and 
subcritical at the individual level with subcriticality parameter eg and with im- 
migration of individuals from Bj^_\ . Moreover by Lemma 15.3.11 the source of 

immigrant individuals Ct+f\t) is asymptotically constant in the time interval 
(t (N),0) provided that 

\to(N)\ ;Q 

as N — > oo. 

The main difference from the proof of Proposition 13.2.11 is that we cannot 
assume that each immigrant belongs to a different family. Since the individuals 
immigrating into Bg N ^ come from families in B^\ subject to family branching, 
it is necessary to keep track of the family structure in the ball B^f[. The reason 
for this is that in principle an individual immigrating into Bg N ^ could be a 
member of a B^_\ -family already represented in B^^ and then could not be 

viewed as the founder of an independent family in B^ N \ Part of the argument 
below is to verify that that this effect is negligible. 

By the family size constraint at time to (see Lemma l5.2.1f b^. the population 
in B^l consists of a collection of families whose sizes are 0{N i - i+r >/ 2 ). We index 
these families at time to by i € N with masses yiN™ +1 >' 2 . 

Recall from section ITO that ty(- N < e )(t) = ^^(f,-) where ipj(t,-) is a 
counting measure on Qn corresponding to the spatial distribution of the family 
indexed by j at time t. We now give a precise formulation to the time develop- 
ment of the families simultaneously in B^_\ and B^ N \ To each family ipj and 
t G K we associate a couple (xj(t),yj(t)) where 

x jVJ- N e/2 ' N (i+i)/2 ■ 



Then 



y {N \t;dx,dy) = e 2 ^5 {Xz(N t,2 t):yr(N i/2 t)) {dx,dy), 



where e 2 = AT-^ +1 )/ 2 . Also, let e = ±, e x = j^. 

By assumption, the family size distribution in B^_\ at time ta, asymptoti- 
cally as N — > oo, is given by fx with J yp?(dy) < oo. 
It suffices to show that for t > to, as N — > oo 

y( N Ht;dx,dy)l {x<N _ i} => n° t+1 (dy)6o(dx) (5.3.5) 

y( N) (t;dx,dy)l {x>0} ^0 (5.3.6) 

{y {N) (t)} t > t0 ^{ V (t)} t > ta (5.3.7) 
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where y( N *> (t) is the renormalized family measure given by 

yW(t;dx,dy) = N^yWfrdx^y^^+yWfrdx^dy)!^.^ 

and where rj(t) is a M 1 [[0, oo) x (0, oo)]-valued diffusion with generator 25^. 

To verify l|5.H.5[l . first using Lemma lf) .2.31 and Proposition ^ . 2 . ll we can verify 
that 



sup 

te(-t ,0] 



(lAy)\y (N) (t;dx,dy)-y W (t o ;dx,dy)\^0 in probability. 



Moreover 

N i/2y(N) (t-dx,dy)<N 1 ' i I I xN l l 2 y( N \t;dx,dy). 



/Af-l/4 

Therefore by the tightness of / / °° x{N 1 / 2 y^ N \-\ dx, dy)l {x>0} ) and {C { e +f } }, 
as N — > oo, 

[ [ y( N \\ Ay) (t;dx,dy) -> in probability 

J Jn- 1 / 4 

and (jfHT5f and (|5~TH| follow. 

We now turn to the proof of (|5.3.7Jl . In order to implement the rescaling we 
introduce the class of functions of the form 

= V?)) = /( / J <P( x i y)v>{dx, dy)) = f(e 2 ^ pfa, Vi)), 



where 



<p(x,y) = e 1/2 l {x>0} ip e (x,y) + ip e+ i{y)l {x<N _i } . 



We assume that \ipe\, | dvi Q^ v ^ \, \ f£+i \ are bounded, C 2 , and \<pi{x, y)\ < const -x 
for x > and </^+i(y)| < const ■ y for y > 0. 
Now define 

u e (dx,dy) = N 1/2 u(dx,dy)l {x>0} 

and note that 

<p(x,y)(i(dx,dy) = J <p e {x,y)ne(dx,dy) + J Lp e+1 (x,y)l {x<N _i } u(dx,dy). 
The generator of y( N ^ acting on F is given by 
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31=0 J 2 >31 

OO OO / • ' \ 
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V(jl£l,j 2 E 2 ) 



2 f-' £ T i . £l£2 
31=0 3 2 >31 

oo oo 

E E 

J1 =1 J2>J1 

JiM(ji£i.J2£2) 
£l£2 

oo oo 
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•[/«/*,¥>> -£l^0l£l,j2£ 2 )+£l^(0l - l)£l,(i2 - 1)£ 2 ))) -/((HM)]^ 161 '^ 



N ~ 



£l£2 
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£l£2 
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\ E E 
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■[f ({H, Vl ) - enptfrexM + e lVt (jiei, (j a - l) Ea ))) - f((n, , 
We give a brief explanation of these terms. 

1. This term corresponds to the migration into b\ N ^ of individuals from 

families currently minimally occupying b\ N \ Asymptotically as N — > oo, this 
is 

ctf'{{n,ip)) / / (—<pi(x,y))yfit + i(dx,dy), 
whose limit (by 15.3.51) ) is 



cef'iifaw)) / (^^( a; >y)lx=o) J/p (dj/). 

2. The second term is similar to the first except that here only families 
having more than minimal mass in -Be appear. In the limit N — > oo we have 
asymptotically 

ctf'((v,tp)) J ytj J^^'^jv^/a M,N(dx,dy) (5.3.8) 

< q|/ ((m, </>})! sup |— <^(x,y)| 

^1/8 /-oo ,7V 1 / 8 

y^(dx, dy) + — m I I ni, N (dx, dy) 

ly>Ni N L / Z J N -l/i J I 

d 

< ci\f ((m>¥>))I su p l^-^O^z/)! 

i yn(dx,dy) + N~ 1/8 / x^^ N (dx,dy)\ 

ly>Ns JO JO ' 

= (1) + 0(^L I ) (5.3.9) 

since J J ;c>0 J x m iN (dx,dy) = 0(1). 

3. and 4. These two terms arise from the family level branching. If 
\tpe(x, y)\ < const • x, then asymptotically they yield 

2^Tj2 /"((M.y)) y y f 2 {x,y)^(dx,dy) 



+/ y y <Pe+i{y)<Pi(x,y)ii(dx,dy) 



5. and 6. This corresponds to the critical birth and death of individuals in 

f'((p,w))(p,G&><p t ), 



b[ N \ The limiting term is 
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(2) 

where Gc e denotes the application of the operator G C( to the x variable. 

7. and 8. These terms correspond to the birth and death of individuals in 
Bg^l in families that do not or sparsely occupy B^ N ^ . Asymptotically we obtain 

1 d 2 

-J «M>¥'))(M ! y-5-o^+i(y)) 



N y 2J wruxniQy 



9. and 10. These terms correspond to the birth and death of individuals 
^e+l m families which 1 
Asymptotically we obtain 



in B^_l in families which have at least a minimal number of members in B^ N ^ 



f 1 d 2 1 

/ ((M, ¥>» / l^Mdx, dy)y—<p e (x, y) = 0(— ^) 

Collecting the limiting terms as N — * oo, we obtain &^F(fi) -> <8 e F(n) 
where 

= f\( f i,<p})( m ,G% > <p l ) + ±f"(((i,<p i })(v i ,$) (5.3.10) 
+f'{(^, l Pe))ce J ^-(x,y)\ x =oy^°(dy). 
We conclude that for any limit point of the probability laws of y^ N 'l x >o, 

F(ri(t)) - f ® e F(ri(s))ds 

Jt 

is a martingale where F € C(M 1 ((0, oo))) is defined by F(/i) = f((fi,(pe) + 
(fi°,ipi+i)) and /, y>£ and ipe+i satisfy the same conditions as above. But this 
coincides with the martingale problem of the two level branching diffusion with 
constant multitype immigration source with immigration rate of type y given 
by y^°{dy) and total immigration rate of J °° y^P(dy) (recall Proposition 13.2. lft 
which is well posed and determines a M 1 (0, oo)-valued diffusion process. This 
completes the proof of the proposition. 



5.4 The limiting multiscale transition function 

The main result of the section is the identification of the limiting multiscale 
structure, which serves to complete the proof of Theorem l2.4.1l 

Theorem 5.4.1 Fix j > 1 and > 0. Conditioned on Cj+f\o) = (f+V, 
and provided that £j+i Cj*+i as ^ °°> then the normalized equilibrium 
masses {Ci^ (Q)}^=j+i in a sequence of nested blocks {-Be }^=i ...j+i 
converge as N — > oo in distribution to the backward Markov chain {d}e=j+i,j....,i 
with transition kernel 



P(Q&A\Q +1 = a) = uf\a,A) 



where is as in Provosition \2.1.H b), and Cj+i = Cj+i 
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Proof. The principal step in proving this is given by Proposition l5.3.2l Then to 
verify that the limiting equilibrium distribution as N — > oo of C^'^ conditioned 
on Ci+i^ — o, is given by IT^ 2 ' (o, dx) we then follow the argument in Corollary 
13.2. 2f b). The convergence of {Ce(fy}i=j+i ,j o to the Markov chain then follows 

(2) 

by recursion and the continuity of the mappings a — > II), (a,dx). 
■ 

Combining Theorem l5.4.1l with the spatial ergodic theorem ( Theorem 15 .1.7(1 
and Proposition 12. l.ll b). we see that the interchange of the limits N — > oo and 
j ' ^ oo is justified and completes the proof of Theorem l2.4.1l 



5.5 The particle level picture 

In the previous section the main result is obtained using the convergence of the 
solutions of the appropriate sequence of martingale problems. Some additional 
understanding of the limiting process can be obtained by examining the particle 
picture in both B^_\ and i?j which we now briefly sketch in an informal 
way. First note that the change in the population in B^f[ due to movement of 
particles between B^+l and its exterior in time scale N e ^ 2 has expected value 
o(jj) and therefore is negligible. 

For the moment we fix to < and consider the contributions to C/ N ' 6 \ty 
coming from immigrants arriving in the interval [toN e / 2 , 0). Next we recall that 
by (|5.2.1|l all but O(^) of the population in B^\ at time is contained in 
families with sizes in (Q,KN t - l+1 ^ 2 ). Therefore we can subdivide the popula- 
tion at time into 0(N^ e+1 ^ 2 ) independent subpopulations of size 0(N^- e+1 ^ 2 ) 
where if necessary we group together smaller families to form subpopulations 
of size 0{N^ +1 ^/ 2 ). (This ensures independent level two branching in B^ for 
clusters coming from distinct families in this subdivision.) Since the families 
undergo critical branching these are the descendants of 0(N 1 ^ 2 ) ancestral fam- 
ilies at time toN e / 2 and in the time interval (toN e / 2 ,0) these undergo binary 
branching with njv(s) branches at time sN^/ 2 . Now consider a single family in 
Bg+l containing 0(N^ e+1 ^ 2 ) individuals. By the analogue of Lemma 15. '6. II the 
normalized population size of this family is constant in the time scale N e / 2 . 

We subdivide the interval [N e / 2 t ,0) into M(€ N) subintervals of equal 
length -Lt N e / 2 . Let t k := and consider an interval (t k+1 N £ / 2 , i fe A^ /2 ). 

Recall that the rate of migration of individuals into B^ from B^\ is ciN~ e / 2 . 
Therefore in the time interval the number of individuals to immigrate into the 
ball B^ from a family of size mN^ e+1 ^ 2 in B^\ is Poisson with mean 

= c i m(t k -t k+1 )NV- 1 » 2 . (5.5.1) 

Now consider the question of which of these have descendants alive in fl| at 
time zero. Let U := min{s : njsr(s) = 2}. In order to determine this recall that 
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from the structure of the genealogy of the critical branching cluster (see e.g. 
|Du| . |F ] ) the random variable U is uniform distributed on [0, 1]. Moreover the 
probability that an initial individual produces a non-empty set of descendants 
in a time interval of length U\tk\N 1 / 2 is asymptotically as N — > oo 



lege 



-c e Ut k 



N e / 2 (1 - e - CtUtk )' 



(see e.g. |AN) fChapt. 3)). Therefore the probability that that any of the 

immigrants arriving from a family of size N<- e +V/ 2 \n[t k+1 N l l 2 ,t k N l l 2 ) 
have descendants alive at time is no larger than 

O r , p -C£U (N)t k 

■P(U<u (N)) 



for any < uq(N) < 1. Again by ( |Du| .[F ) ') the number of (family) branches 
^Af(ife) at time t k N l l 2 converges in distribution as N — > oo for each of the 
0{N X / 2 ) ancestral families and the family trees from the different ancestral 
families are independent. Choosing uq(N) — > such that N 1 / 2 uq(N) — > oo 
we conclude that there is a Poisson [with mean of order 0((t k — t k+1 )e~ citk )\ 
number of individuals who immigrate in (N^^to, N e / 2 tk] producing descendants 
at time 0, and these all come from different branching trees. 

Therefore the population at time zero coming from immigrants arriving in 
(N i / 2 ta,N i l 2 tk) is asymptotically composed of 0(1) two level families each 
originating from one individual and these all come from different independent 
subpopulations in B^f^. In particular the expected mass coming from clusters 
containing two or more immigrants from the same subpopulation in B^?\ is 
0{ N \ /2 ). Each two level cluster coming from an immigrating individual develops 
by two level branching, namely, family level branching inherited from the family 
branching in B^+l and subcritical individual level branching in B^ N \ Finally, 
in both cases the total migration rate of individuals into B^^ from B^f\ is 

(C!+i e) (A^o)iV £+1 ) x (c e N-^ 2 ) x TV- 1 = c e N^ 2 C^f\N^ 2 t ) 

where on the left hand side the first factor is the number of particles in B^_\ , 

the second is the individual migration rate to a point chosen randomly in B^[\ 
and the last factor is the probability that the tagged ball is chosen. There- 
fore, asymptotically as N — ► oo, the population at time zero consists of clusters 
of descendants of individuals that immigrate into B^ during the time inter- 
val (—oo,0] at rate C£N e ^ 2 (^f\N e ^ 2 t ) and subsequently undergo two level 
branching. In the limit N — > oo these clusters correspond to the jumps of the 
subordinator Si(-) defined in section POl 

Remark 5.5.1 One can also gain some understanding of the convergence to 
equilibrium from a spatially homogeneous initial population (more general than 
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that addressed in Proposition \2.S.S6 with intensity 9 > 0. Two ingredients 
are involved in the convergence to an equilibrium with intensity 6. The first 
is the strong transience condition on the random walk. The other feature is 
the structure of the local family sizes. We see from the above that the property 
that the contribution of families in B^J containing a number of individuals 
larger than 0(N^ +1 ^ 2 ) to the equilibrium population in B^^ is asymptotically 
negligible and this property is then inherited by B. . However if the initial 
family sizes are too large this iteration can degenerate due to the family level 
critical branching and the limiting population is locally degenerate. For a more 
detailed analysis of this phenomenon, see ]GHy . 



6 Appendix 

6.1 Size-biasing and Palm distributions 

Definition 6.1.1 If n is a measure on some measurable space M, and s is a 
nonnegative measurable function on M with < s := J s{z)Tt(dz) < oo, then we 
call the probability measure tt given by 

Tt(dz) := ^ir(dz) 
s 

the size-biasing of it with respect to 5. (Here, we think of s(z) as measuring the 
size of the object z.) 

Remark 6.1.2 a) An example of size-biasing which is important in our context 
arises as follows. Let n — C(rj) be the distribution of a random measure n (on 
some Polish space E, say). Denote the intensity measure of r) by X, and fix a 
nonnegative measurable function f on E with < (A,/) < oo. Define the size 
of a measure m on E by s(m) := J f(x)m(dx) — (m, /}, and denote by itf the 
size-biasing of it(dm) with (m, f) . 

b) Assuming that A is locally finite, choosing f = 1b where B is a ball in E, 
and letting B shrink gives the family of Palm distributions n x , i£5. Formally, 
these arise as the disintegration of the measure n(dm)m(dx) with respect to its 
second marginal E?y 7 that is 

EG(n){r),h) = J h(x)G(m)n x (dm)\(dx). (6.1.1) 

See fffi . chapter 10, for more background; there, a random measure whose dis- 
tribution is the size-biasing of (C(rj))(dm) with (m, /) is denoted by r)f. 

c) The following fact (JK^. formula (10.6)) is immediate from jd.l.lp : 
The size-biasing TTf of it (dm) with (m, /) is 

f(x)ir x (.)\(dx) =E7r jt , (6.1.2) 



(A,/) 
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where X is a random element in E whose distribution is the size-biasing of A 
with f. 

d) If in the just described situation E consists of one element only, then the 
finite random measures on E are R + -valued random variables. When speaking 
of the size-biasing of a measure -k on R + without specifying the size function, 
we always mean the size-biasing of ir(dx) with x. 

Let us write II a for the distribution of a Poisson random counting measure on 
E with intensity measure A. It is well known (see e.g. jK], beginning of chapter 
11) that the Palm distributions of 11a arise as the distributions of Q + 8 X , x G E, 
where £($) = IIa. 

Now let a be a probability measure on R + with ra a :— J Ta(dr) G (0, oc), 
and write 

IT^a := J U T \a(dT) 

for the mixed Poisson distribution with mixing measure a. 

The following lemma, whose proof we include for convenience, is part of a 
characterization theorem ([K], Theorem 11.5) of mixed Poisson processes. 

Lemma 6.1.3 The Palm distributions of Ilcr, a arise as the distributions o/$ + 
8 X , x G E, where C(Q>) = Tla.x, and a is the size-biasing of a. 

Proof. Using the above mentioned form of the Palm distributions of II t a, we 
obtain for all nonnegative measurable h and G defined on E and M C (E), the 
space of locally finite counting measures on E, respectively: 

G(^)^ft)n M (#) = [ [ G(iP)(iP,h)Il TX {di>)o-(dT) (6.1.3) 

-\{dx)h{x)G{ip + 5 x )U TX {d^)a{dT) 

m a \(dx)h{x) ( G(tp + S x )H T x(dtp) raidr) 



Corollary 6.1.4 Assume < (A, h) < oo. Then the size-biasing of Tl at \(dil)) 
with (ip,h) arises as the distribution of&+6j{, where £(<&) = J H T \a(dT), C(X) 
is the size-biasing of X with h, and $ and X are independent. 

Proof. This can be seen either by combining Remark lG. 1.21 c) and Lemma r6.1.3l 
or by dividing (|6.1.3|l through E CTi a( v 1', h) = m CT (A, h). ■ 



6.2 Subcritical Feller branching 

Let us fix c > 0. In the following, X will denote a c-subcritical Feller branch- 
ing diffusion (c-FBD) process. In other words, X is an [0, oo)-valued diffusion 

process satisfying 

dX t = x/YtdWt - cX t dt (6.2.1) 
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where Wt is Brownian motion. For e > 0, let X s be the c-FBD process starting 
in e at time 0. 

From well-known results on Galton- Watson processes conditioned to sur- 
vival jANHGej and cluster sizes in continuous-state branching processes [D], one 
expects that the conditional law C(Xf\Xf > 0) converges to an exponential 
distribution as e — > 0. The following calculation verifies this and identifies the 
parameter. 

The Laplace functional of Xf is given by 

E( e - AX *) = e - ev{t - x \ A>0, (6.2.2) 
where v = v(t, X) is the solution of 

^L^l^^cv(t,X)-^v 2 (t,X), v(Q,X)=X. (6.2.3) 
The solution of (|6.2.3[1 is given by 

2Ace 
A(l -e- ct ) + 2c 

Combining (|6.2.2|l and l|6.2.4(l one obtains by a straightforward calculation 



v(t,X)= W1 __ et . , - ■ (6.2.4) 



F(Xf ± 0) = 1 - exp y x _\_ ct j , (6.2.5) 

EX t e = ee~ ct , e > (6.2.6) 

and 

lim E(e- AX * \Xf ± 0) = — ^ — . (6.2.7) 

e^o v A(l - e- c *) + 2c v 7 

Writing Exp(u) for the exponential distribution with parameter u, we obtain 
immediately: 

Lemma 6.2.1 Fix t > 0. 

a) 

2c 
1 ~~ 



L(X e t \X e t > 0) => Exp ( — ) (6.2.8) 



2ce" ct 



1 -e 
For e > 0, let 



X P[X^ > 0] -> — as e -> 0. (6.2.9) 



X e = ^X £ <\ (6.2.10) 



where 7V e is a Poisson(e 1 )-random variable and X £,1 : X 6,2 , . . . are independent 
copies of X e . The following lemma is an easy consequence of Lemma 16.2.11 
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Lemma 6.2.2 Fix t > 0. Then 

X e t X° as e -+ 0, (6.2.11) 
where X® is infinitely divisible with canonical measure 

Let us now define the measure 

j c := c / K s ds. (6.2.13) 
Jo 

An elementary calculation based on H6.2.12J) and H6.2.13J) shows that 

lc {dx) = 2c-e~ 2cx dx, x>0. (6.2.14) 

x 

This identifies j c as the canonical measure of the Gamma(2c, 2c)-distribution, 
and goes along with the well-known fact that the equilibrium distribution of 
(|4. 1.1(1 is the Gamma(2ca, 2c)-distribution. We note that 7 C is the Levy measure 
of the Gamma subordinator S(t), t > 0, with scale parameter 2c and ES'(l) = 1. 
The following is obvious from l|6.2.14|) : 

Remark 6.2.3 a) f °° xj c (dx) = 1. 
b) f™ x 2 j c (dx) = l/2c. 

Finally we consider the semigroup (T t ) be the semigroup of the c-FBD pro- 
cess. Recalling 1)6.2.80 . (|6.2.9[) and l|6.2.12jl . the (T t )-entrance law (re t ) from is 
given by 

K t = lim -8 e T t = lim -£(Xf; Xf ^ 0), t> 0, (6.2.15) 
e—>0 e e— >o e 

where (Xf) is the c-FBD-process starting in e at time 0. Because of <|6.2.12ll . 
K t (dx) has density 

(2c) 2 e~ ct / 2cx 

Then for / 6 Ci((0, oo)) and t > 

d_ 

dx 



K t( x ) = f) \-cn2 ex P ( - 1 — — ) . 1 e (°. ( 6 - 2 - 16 ) 



T t /(x) =lim-T t f(e)= [ f(y)K t (dy). (6.2.17) 
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